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Abstract. Our main result is a limit shape theorem for the two-dimensional surface 
defined by a uniform random n x n square Young tableau. The analysis leads to a calculus 
of variations minimization problem that resembles the minimization problems studied 
by Logan-Shepp, Vershik-Kerov, and Cohn-Larsen-Propp. Our solution involves methods 
from the theory of singular integral equations, and sheds light on the somewhat mysterious 
derivations in these works. An extension to rectangular diagrams, using the same ideas 
but involving some nontrivial computations, is also given. 

We give several applications of the main result. First, we show that the location of a 
particular entry in the tableau is in the limit governed by a semicircle distribution. 

Next, we derive a result on the length of the longest increasing subsequence in seg- 
ments of a minimal Erdds-Szekeres permutation, namely a permutation of the numbers 
1, 2, . . . , n 2 whose longest monotone subsequence is of length n (and hence minimal by the 
Erdos-Szekeres theorem). 

Finally, we prove a limit shape theorem for the surface defined by a random plane 
partition of a very large integer over a large square (and more generally rectangular) 
diagram. 
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1 Introduction 



1.1 Random square Young tableaux 

In this paper, we study the large-scale asymptotic behavior of uniform random Young 
tableaux chosen from the set of tableaux of square shape. Recall that a Young diagram 
is a graphical representation of a partition A : A(l) > A(2) > . . . > X(k) of n = J2 ^ as 
an array of cells, where row i has A^ cells. For a Young diagram A (we will often identify 
a partition with its Young diagram), a Young tableau of shape A is a filling of the cells 
of A with the numbers 1, 2, . . . , n such that the numbers along every row and column are 
increasing. 

A square Young tableau is a Young tableau whose shape is an n x n square Young 
diagram. The number of such tableaux is known by the hook formula of Frame-Thrall- 
Robinson (see © below) to be 

(rf)\ 

[1 • (2n - 1)][2 • (2n - 2)] 2 [3 • (2n - 3)] 3 ...[(n- l)(n + l)]"- 1 n n ' 

A square tableau T — (tij)2j = i can be depicted geometrically as a three-dimensional 
stack of cubes over the two-dimensional square [0, n] x [0, n], where cubes are stacked 
over the square [i — X \j — 1, j] X {0}. Alternatively, the function — » fcy can 
be thought of as the graph of the (non-continuous) surface of the upper envelope of this 
stack. By rescaling the n x n square to a square of unit sides, and rescaling the heights of 
the columns of cubes so that they are all between and 1, one may consider the family 
of square tableaux as n — > oo. This raises the natural question, whether the shape of the 
stack for a random n x n square tableau exhibits some asymptotic behavior as in oo. 
The answer is given by the following theorem, and is illustrated in Figure 1. 



Theorem 1. Let T n be the set ofnxn square Young tableaux, and let P„ be the 
uniform probability measure on T n . Then for the function L : [0, 1] x [0, 1] — > [0, 1] defined 
below, we have: 

(i) Uniform convergence to the limit shape: for all e > 0, 



\[TE%: max 

l<i,j<n 



f 



* J 



> e 



0. 



(ii) Rate of convergence in the interior of the square: for all e > 0, 



TeT„: max 

1 < i, j < n 
min(ij, (n — i){n — j)) > n"^ 2- *" 6 
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(a) 3D plot of simulated tableau (b) The limit surface L(x, y) 




10 20 30 40 50 0.2 0.4 0.6 0.8 1 

(c) Contour plot of simulated tableau (d) Contour plot of L 

Figure 1: A simulated 50 x 50 random tableau and the limit surface 
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Definit ion of L. We call the function L the limit surface of square Young tableaux. 
It is defined by the implicit equation 

2 _ x ( (l-2L{x,y))(x-y) \ 

x + y = — (x — y) tan 1 1 



7T- -■ \y/4L(x,y)(l-L(x,y)) - (x - y) 2 J 

2 x / ^/4L(x,y)(l-L(x,y))-(x-y)A 
tt an ^ l-2L(x,y) J 

for < y < 1 — x < 1, together with the reflection property 

L(x,y) = 1 - L(l - x, 1 - y) 

(where tan -1 is the arctangent function). It is more natural to describe L in terms of its 
level curves {L(x,y) — a}. First, introduce the rotated coordinate system 

X-y x+y 

In the u — v plane, the square [0,1] x [0, 1] transforms into the rotated square 

O = {(u, v) e R 2 : M < y/2/2, \u\ < v < V2 - |u|}. 
Now define the one-parameter family of functions (g a )o<a<i given by 



g a : [-y/2a(l-a), y/2a(l - a)] 



g a {u) 



l^tan- 1 ( , {1 - 2a ^ ) + ^tan' 1 ( 0< a <± 

* \ V2a(l-a)-u 2 J » \ 1 ~ 2a J - 2 

-lutan-i ( , {1 ~ 2a)u , ) - ^tan^ ( V^ad-a)- ^ \ ^ i (2) 

"~ \ v / 2a(l-a)-u 2 y 71 \ 1 ~ 2a J ~ 

2 " 2 ■ 



Then in the rotated coordinate system, the surface L{u,v) = L(x(u,v),y(u,v)) can be 
described as the surface whose level curves {L(u, v) — a} are exactly the curves {v = 
g a (u)}. That is, 



{(u,v) € O : L(u,v) = a} = {(u,v) £ O : \u\ < y/2a(l - aj,v = g a (u)}. 

This is illustrated in Figure 2. It is straightforward to check that the curves v = g a (u) do 
not intersect, and so define a surface 2 . 



2 See equation (ff)f))> in section 3.4. 
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Figure 2: The curves v = g a (u) for a = 0.05, 0.1, 0.15, 0.2, . . . , 0.5 



Note some special values of L(x,y) which can be computed explicitly: 



L(t,0) = L(0,t) = 1 V T~ f2 



r , \ w- x 1 + \f2t - t 2 
L(t,l) = L(l,t) = y — . 

= l-COsCTrf)^ 



The approach in proving Theorem 1 is the variational approach. Namely, we identify 
the large-deviation rate functional of the level curves of the random surface defined by 
the tableau, then analyze the functional and find its minimizcrs. This will give Theorem 
l(ii), with the rate of convergence following from classical norm estimates for some integral 
operators. The treatment of the boundary of the square, required for Theorem l(i), turns 
out to be more delicate, and will require special arguments. 



1.2 Location of particular entries 

Theorem 1 identifies the approximate value of the entry of a typical square tableau in a 
given location in the square. A dual outlook is to ask where a given value k will appear 
in the square tableau, since all the values between 1 and n 2 appear exactly once. These 
questions are almost equivalent. Indeed, if k is approximately a ■ n 2 , then Theorem 1 
predicts that with high probability the entry k will appear in the vicinity of the level curve 
{L(x, y) = a} (the fact that this actually follows from Theorem 1 is a simple consequence 
of the monotonicity property of the tableau along rows and columns). However, one may 
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ask a more detailed question about the limiting distribution of the location of the entry 
k on the level curve. It turns out that its it-coordinate has approximately the semicircle 
distribution. This is made precise in the following theorem. 



Theorem 2. For a tableau T e T n and 1 < k < n 2 , denote by [i(T,k),j(T,k)) 
the location of the entry k in T, and denote X(T,k) = i(T,k)/n, Y(T,k) = j(T,k)/n. 

Let < a < 1, let k n be a sequence of integers such that k n /n 2 > a, and for 

each n let T n be a uniform random tableau in T n . Then as n — > oo, the random vector 
(X(T n , k n ), Y(T n , k n )) converges in distribution to the random vector 

(X a ,Y a ) :- ^ , — 

where U a is a random variable with density function 



, , x ^2a(l-a)-u 2 

/l/ ° (u)= ^(1 -a) VVmI^V^M^] {u) (3) 

and V a = g a (U a ). 

Theorem 2 is one of several aspects of our work which shows a deep connection to 
the work of Logan-Shepp and Vershik-Kerov on the limit shape of Plancherel-random 
partitions - see section 8 for discussion. 



1.3 Minimal Erdos-Szekeres permutations 

The famous Erdos-Szekeres theorem states that a permutation of 1,2, ... ,n 2 must have 
either an increasing subsequence of length n or a decreasing subsequence of length n. This 
can be proved using the pigeon-hole principle, but also follows from the RSK correspon- 
dence using the observation that a Young diagram of area n 2 must have either width or 
height at least n. 

For the width and height of a Young diagram of area n 2 to be exactly n, the diagram 
must be a square. From the RSK correspondence it thus follows that to each permutation 
of 1,2, ... ,n 2 whose longest increasing subsequence and longest decreasing subsequence 
have length exactly n, there correspond a pair of square n x n Young tableaux. Such a 
permutation has the minimal possible length of a longest monotone subsequence, and it 
seems appropriate to term such permutations minimal Erdos-Szekeres permutations (we 
are not aware of any previous references to these permutations, aside from a brief mention 
in exercise 5.1.4.9). 

As an application of our limit shape result, we will prove the following result on the 
length of the longest increasing subsequence when just an initial segment of a random 
minimal Erdos-Szekeres permutation is read. 
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Theorem 3. For each n, let w n be a uniform random minimal Erdos-Szekeres per- 
mutation of 1, 2, ... , n 2 . For 1 < k < n 2 , Let be the length of the longest increasing 
subsequence in the sequence 7r„(l), 7r„(2), . . . , 7r„(fc). Denote a — k/n 2 , and cto = n~ 2 / 3 + c . 
Then for any e > 0, and w(n) — * oo however slowly, 

max V(\l n k — 2-\/q!(1 — a)n\ > a^ 2 uj(n)n) > 0. 

a Q <k/n 2 <l/2 ' n— «x> 

Thus the random fluctuations of l n .k around 2^Ja{l — a)n are not likely to be of order 
substantially larger than n 2 / 3, . 

1.4 A limit surface for random square plane partitions 

Another probability model which was studied in the context of limit shapes, is that of 
random plane partitions. If A is a Young diagram, a plane partition of n of shape A is an 
array of positive integers {pij)(ij)£\ indexed by the cells of A, that sum to n and which 
are weakly decreasing along rows and columns, i.e., satisfy 

Pi,j > Pi,j • i- Pid > Pi+i.j- 

Cerf and Kenyon [7j proved a limit shape result for random unrestricted plane partitions 
of an integer n. Cohn, Larsen and Propp 8 proved a limit shape result for random plane 
partitions whose three-dimensional graph is bounded inside a large box of given relative 
proportions. 

We apply Theorem 1 to prove a limit shape result for random plane partitions of m 
defined over an n x n square Young diagram, when m is much greater than n . This can 
be related to square Young tableaux by the observation that when a plane partition does 
not contain repeated entries (which in the asymptotic regime described above happens 
with high probability), the order structure on the entries of the plane partition is a Young 
tableau. The precise result is the following. 

Theorem 4. For integers n, m > 0, let V n ,m be the set of plane partitions of m of n x n 
square shape, and let P n ,m be the uniform probability measure on V n , m - If 7r = (j>i,j)i l j = i is 
an element of V n ,mi let its reseated surface graph be the function : [0, 1) x [0, 1) — > [0, oo) 
defined by 

^T:(x,y) = —P[nx]+l,lnyi+l- 

Suppose m and n are sequences of integers that tend to infinity in such a way that 
m/n 6 — * oo. Then for all e > 0, x, y £ [0, 1) we have 

y)-M(x,y)\ > e) ► 0, 
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where M : [0, 1] x [0, 1] — > [0, oo) is given by 

M{x,y) = -log(L(x,|/)). 

Theorem 4 may be related to a limiting case 7 — ► 00 in the limit shape result of 
Cohn-Larsen-Propp |SJ. We have not attempted to check this. 

1.5 Random rectangular Young tableaux and plane 
partitions 

The methods which we will use to prove Theorems 1, 2, and 4 work equally well for 
rectangular Young tableaux and plane partitions, in the limit when the size of the rectangle 
grows and its relative proportions tend to a limiting value 9 > 0. For each possible value 
9 of the ratio between the sides of the rectangle, there is a limiting surface Lg for random 
rectangular Young tableaux, and a limiting surface Mg for random rectangular plane 
partitions. Analogously to the square tableaux, the rectangular m x ri2 tableaux can be 
viewed as the result of applying the RSK algorithm to a permutation of {1, . . . , n\n2\ 
with the property that the lengths of the longest increasing and the longest decreasing 
subsequences are exactly equal n\ and ni (by the Erdos-Szekeres theorem, the two lengths 
cannot be simultaneously below n\ and respectively). The proofs, which we include 
at the end of the paper, require some nontrivial modifications, but the final results are 
unexpectedly as elegant as for the square case. 

Let 9 > 0. We may assume that 9 < 1, otherwise exchange the two sides of the 
rectangle. Define Lg : [0, 1] x [0, 9} — > [0, 1], the limit surface of rectangular tableaux with 
side ratio 9, as follows. For each < a < 1, the a-level curve {(x,y) : Lg(x,y) = a} is 
given in rotated u — v coordinates by 

{(u, hg lCt (u)) : -f3i <u< f3 2 }, 
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where 



(3 = ^28a{\ - a), 

(3i = 0-a(l -0)V2/2, fa = P + a(l-9)V2/2, 



h e , a {u) = 9V2/2 ± (/?! - 9V2/2) + 



± (-£-7i) tan 



(1-0(71-1) 
V (1 + C)(7i + 1) 



+ (C-72)tan" 



.! (1 + 0(72-1) 



1 / • 1 - T 



(1-0(72 + 1) 

7T\ 1 -g 7T ' 

Vf>f ± 2 (7l - 1) 



0< a< ^> 



± = 



7i 

he, a {u) 
see Figure 3. Set 



J + < a < 0/(1 + 0), 
\ - 0/(1 + 6) <a<l/2 

u - a(l - 6>)V2/2 



/3 

a + 6»(1 - a) 



(1 + 6<)V2/2 - /i M -a((l - ^)V2/2 - u), 
Me (a;, y) = - log(L e (a:, j/)). 



8a + 1 — a 



2<«<1. 



Theorem 5. For integers n,m > 0, let 7^ jm be the set of tableaux whose shape 
is an n x m rectangular diagram, and let W n ,m be the uniform probability measure on 
P n ,m- If T = (Uj)ij G T n ^ m , define the rescaled tableau surface of T as the function 
S T : [0, 1) x [0, m/n) -> [0, 1] given by 



Sr{x,y) 



1 



nm 



*|na:J+l,|n2/J+l- 



Let < 6 < 1. If m„ is a sequence of integers such that m n /n — > as n — > 00, then for 
all e > 0, x e [0,1), y e [0,9), 



11 n,m n (T E T n , mn : |5 T (x,y) - L e (x,y)| > e) 



0. 



Theorem 6. For each 8 > 0, Mg is the limit surface of uniform random plane partitions 
of m„ over a rectangular diagram of sides n and k n , provided k n /n — > and m n /n 6 — > 00 
as n — > 00. (The precise statement is by analogy with Theorems 1, 4, and 5.) 
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Figure 3: The curves he, a for = 0.5, a = k/9, k = 1, 2, . . . , 8. 
1.6 Organization of the paper 

The remainder of the paper is organized as follows: In the next section, we present the 
variational approach to the limit surface of random square Young tableaux, based on the 
hook formula of Frame-Thrall-Robinson. The level curves of L appear as minimizers of a 
certain functional. This leads to a proof of Theorem 1 in the interior of the square, except 
for the explicit identification of L. Section 3 is dedicated to the derivation of the explicit 
formula for the minimizer. Unlike the analogous results of Logan-Shepp, Vershik-Kerov 
and Cohn-Larsen-Propp, we will show that there is no need to guess the minimizer, by 
giving a technique for systematically deriving it using an inversion formula for Hilbert 
transforms on a finite interval. This may prove useful in similar problems. 

In section 4, we complete the proof of Theorem 1, treating the more delicate case 
of the boundary of the square, and prove Theorem 3. In section 5, we discuss the hook 
walk of Greene-Nijenhuis-Wilf and the concept of the co-transition measure of a Young 
diagram. Using the explicit formulas for the co-transition measure derived in |3(J|. we 
compute the co-transition measure of the level curves g a , proving Theorem 2. In section 
6, we prove Theorem 4. In section 7 we give the computations necessary for settling the 
rectangular case. In section 8, we discuss the connections of our results to the theory of 



11 



Plancherel-random partitions, and some open problems. 

2 A variational problem for random square 
tableaux 

2.1 A large-deviation principle 

One may consider a tableau T G T n as a path in the Young graph of all Young diagrams, 
starting with the empty diagram, and leading up to the n x n square diagram, where 
each step is of adding one box to the diagram. Identify T with this sequence A^ — ^ c 
Ay C C ... C A^ 2 = □„ of diagrams. (A^ is simply the sub-diagram of the square 
comprised of those boxes where the value of the entry of T is < k.) Theorem 1 is then 
roughly equivalent, in a sense that will be made precise later, to the statement that for 
each 1 < k < n 2 — 1, the rescaled shape of A^ for a random TeT„ resembles the sub-level 
set 

{(x,y)e[0,l] 2 :L(x lV )<k/n 2 } 

of L, with probability 1 — o(l) as n — > oo. It is this approach that leads to the large- 
deviation principle. Namely, we can estimate the probability that the sub-diagram A^ has 
a given shape: 

Lemma 1. For T e T n , denote as before A^ C . . . C \^ the path in the Young 
graph defined by T, and for each < k < n 2 , let : A£(l) > A^(2) > . . . > A^(n) 
be the lengths of the columns of A^ (some of them may be 0). For any Young diagram 
A : A(l) > A(2) > . . . > X(n) whose graph lies within the n x n square, define the function 
h ■ [0, 1] - [0, 1] by 

h(x) = h(\nx]). (4) 

(Note that this depends implicitly on n.) Let < k < n 2 , and let a — k/n 2 . Then for any 
given diagram Ao C □„ with area k, we have 

P„ {TeT n :\ k T = A ) = exp ( - (1 + o(l))n 2 (/(/ Ao ) + H(a) + C)) (5) 
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as n — ► oo, where 

C = --21og2, 
H(a) = —a log(a) — (1 — a) log(l — a), 

1(g) = log\g(x) - y + g~ l {y) - x\dydx, 

Jo Jo 

g-\y) = inf{x e [0, 1] : g(x) < y}. 
The o(l) is uniform over all Ao and all < k < n 2 . 

Proof. For a Young diagram A : A(l) > A(2) > . . . > A(7) of area m, denote by d(X) the 
number of Young tableaux of shape A (also known as the dimension of A, as it is known 
to be equal to the dimension of a certain irreducible representation corresponding to A of 
the symmetric group of order m). Recall the hook formula of Frame-Thrall- Robinson 
which says that d(X) is given by 

d(A) = n ^ h ' (6) 

where the product is over all boxes (i,j) in the diagram, and hi a is the hook number of a 
box, given by 

= 1 + number of boxes either to the right of, or below (i,j) 



(and where A' is the conjugate partition to A.) Then we have 



3 



F n (Ter n :4 = X ) = ^^^ , (7) 

where c?(D„\ Ao) means the number of fillings of the numbers 1, . . . , n 2 — k in the cells of the 
skew- Young diagram □„ \ Ao that are monotonically decreasing along rows and columns. 
This is because □„ \ Ao can be thought of as an ordinary diagram, when viewed from the 
opposite corner of the square. The number of square tableaux whose fc-th subtableau has 
shape Ao is simply the number of tableaux of shape Ao, times the number of fillings of the 
numbers k + 1, k + 2, . . . , n 2 in the cells of □„ \ Ao that are monotonically increasing along 
rows and columns - and these are of course isomorphic to tableaux of shape □„ \ Ao , by 
replacing each entry i with n 2 + 1 — i. 



J Note to the reader: this is probably the most important formula in the paper! 
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Take minus the logarithm of |J7J| and divide by n 2 , using ©. The right-hand side 
becomes 

I / ( 2 \\ \ 1 " _ _ 
a + 6 + c-d:=^Iog( ij ^^) + ^X;X;iog(A(0-i + V0Vi + l) 



E log(j-A(i)+i-A'(i) + l)--2^^1og(2n-i-i + l). 



i=l j=X(i) + l i=l j=l 



By Stirling's formula, we have a = nT 2 log 0\ ) = H (a)+o(l), with the required uniformity 
in fc. The other summands look like Riemann sums of double integrals. Indeed, we claim 
that 

1 /-/xoW / _ \ k 

log ( f Xo (x) -y + fx (y) -xjdydx + logn + o(l), 

/ \ n 2 ~ k 

log y - Ao (x)+x- f7 (y) \dydx + — log n + o(l ) , 

/o J/a (x) V / n 

i-l rl 

log (2 — a; — da; + log n + o(l) = C + log n + o(l), 



/o Jo 

»1 /•! 



which on summing and exponentiating would give the lemma. Let us prove, for example, 
the first of these equations. Write 



1 n A« 



= 1 3 = 1 



1 v-v^i / A(z)-i + A^)-z + l \ fc 
= ~Z^Z^ lo M +— logn. 

i=l 3=1 V / 

Fix 1 < i < n and 1 < j < A(i). Denote /i = (A(i) — j + A'(j') — i + l)/n. Approximate 
n~ 2 log h in the above sum by the double integral 

Q ■= / / log ( Ao (x)-y + /r 1 (y) - x) dx. 

J(i-l)/nJ(j-l)/n 

A change of variables transforms this (check the definition of Ao) i n t° 

/ l/2ra /-l/2n 
/ log(a; + y + h)dx dy. 
-l/2n J-l/2n 

Note that ft may take the values 1/n, 2/n, . . . , (2n — 1) /n. If h = 1/n, then integrating we 
get 



»i /.i 

Q = -^f + n- 2 I / log(u + i;)duefo = + 0(n 
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If h > 2/n, by the integral mean value theorem, we have for some r\ G [—1,1], 
Q= l g(/ t + ^) = log^ + 

Clearly then the last estimate holds for h = 1/n as well. The sum of the remainders over 
all 1 < i < n, 1 < j < X(i) is of order 



n 



2n-l , > 

2 y J_ <n - 2 y^nl_ 

ha ~ ^ m 

(i,j)e\ lJ m=l 



where 

a(m) := € A : = m}. 

Clearly a(m) < n, since each row i of Ao contains at most one cell (i, j) with hij = m. 
This gives that the sum of the remainders is of order 

2n-l /, 
2 \ ^ n r\( i°S n 



^ — ' m 

which is indeed o(l). 



n 

m 

m— 1 



2.2 Two formulations of the variational problem 

Lemma 1 says, roughly, that the exponential order of the probability that a random square 
tableau T has a given fc-subtableau shape, where k is approximately a ■ n 2 , is given by the 
value of the functional I on the boundary g of the shape, plus some terms depending only 
on a. Following the well-known methodology of large deviation theory, the natural next 
step is to identify the global minimum of / over the appropriate class of functions, or in 
other words to find the most likely shape for the adevel set. If we can prove that there is 
a unique minimum, and identify it, that will be a major step towards proving Theorem 1. 
So we have arrived at the following variational problem. 

Variational problem 1. For each < a < 1, any weakly decreasing function 
/ : [0, 1] -» [0, 1] such that J* Q f(x)d x = a is called a- admissible. Find the unique a- 
admissible function that minimizes the functional 

/(/)= / / log \f(x)-y + f-\y)-x\dydx. 
Jo Jo 

We now simplify the form of the functional /, by first rotating the coordinate axes by 
45 degrees, and then reparametrizing the square by the "hook coordinates" - an idea used 
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in |201- Let u, v be the rotated coordinates as in QJ. Given an a-admissible 

function / : [0, 1] — > [0, 1], there corresponds to it a function g : \—«j2/2, y/2/2] — > [0, \/2], 
such that 

y = /O) « = ff(w) 

(see Figure 4). The class of a-admissible functions translates to those functions g : 
[-• v/2/2, -s/2/2] -> [0, V2] that are 1-Lipschitz, and satisfy g(-y/2/2) = g(\/2/2) = V2/2 
and 



(<?(it) — |u|)t2tt = a. 



(8) 



'-V5/2 

We continue to call such functions a-admissible. We call a function admissible if it is 
a-admissible for some < a < 1. 




-V2/2 



s t y/2/2 



Figure 4: The rotated graph and the hook coordinates s, t 



To derive the new form of the functional, write 

1(f) = h(f) + Hf) := f [ n \og(h f (x,y))dydx+[ f log(h f (x,y))dydx, 

J0 JO Jo J f{x) 

where hf(x,y) is the hook function of /, 

h f (x ) y) = \f(x)-y + f- 1 (y)-x\. 

Now, set 

J(g) = Ji(g) + Mq) ■= h(f) + W), 
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loghf(x, y)dv du. 



where / and g are rotated versions of the same graph as in Figure 4. Then 

/■a/2/2 pV2-\u\ 
I-V2/2 Jg(u) 

Reparametrize this double integral by the hook coordinates s and t, 

_ r 1 {y) - y , _ x- f(x) 

(see Figure 4). The Lipschitz property ensures that this transformation is one-to-one from 
the region 

{(u,v) : -V2/2 <u< V2/2, g(u) <v<V2- \u\} 

onto the region 

A = {(s, t) : -V2/2 < s < t < V2/2}. 
Therefore the integral transforms as 

<9(it, v) 



Mf) = J J log(V2(*-«) 



d( s ,t) 



ds dt. 



It remains to compute the Jacobian d(u,v)/d(s,t). An easy computation gives (see |35| . 

d(u, v) 1 . „ ,., .... 

(This can be viewed geometrically as follows: draw on the u-axis in Figure 4 the two 
intervals [s, s + ds], [t, t + dt]. The set of points in the square for which the hook coordinates 
fall inside the two intervals is approximately a parallelogram whose area is clearly seen 
from the picture to be linear in 1 — g'(s) and in 1 + g'(t).) So 

' ' r log(V2(t-s)) (l-g'(s))(l+g'(t))dsdt. 



A similar computation for J\, using "lower" instead of "upper" hook coordinates, shows 
that 



Ji(ff) = - / / log(V2(i-a)) (l+g'(s))(l-g'(t))dsdt. 



This gives 

J(g) = 



log (V2(t - a)) [(1 - g'(s))(l + g'{t)) + (1 + g'(s))(l - g'(t))] dsdt 
log U/2(t - sfj (2 - 2g'{s)g'(t))dsdt 

x r V2/2 ,V2/2 3 

-/ / log\t-s\-g'(s)g'(t)dsdt + \og2--. 
z J-V2/2 J-V2/2 1 



We can now state a reformulation of the original variational problem. 
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Variational problem 2. For each < a < 1, a function g : [— v/2/2, V2/2] — * 
[0,a/2] is called a-admissible if: g(-y/2/2) = g(V2/2) = y/2/2 ; 5 is 1-Lipschitz; and 
f_^/ 2 (9{ u ) ~ \ u \)du — a. Find the unique a-admissible function that minimizes the 
functional 

K{g) = ~- / g'{s)g'{t)\og\s-t\dsdt. (9) 

z JV2/2 J-t/2/2 

2.3 Deduction of Theorem l(ii) 

In the next section, we prove the following theorem. 

Theorem 7. For each < a < 1, let cji Q be the unique extension of <7 a (defined in 
to an a-admissible function, namely 



9a(u) 

for < a < 1/2, and 

9a{u) = 



g a (u) \u \ < yggl - q) 

|m| ^/2a(l - a) < |«| < V2/2 



g a (u) \u\ < y/2a(l - a) 
V2-\u\ ^2a{l - a) < |u| < y/2/2 



for 1/2 < a < 1. Then: 

(i) g a is the unique solution to Variational problem 2; 

(ii) K{g a ) = -H(a)+\og2; 

(iii) For any a-admissible function g we have 

>K(g a )+K(g-g a ). 



Assuming this as proven, our goal is now to prove Theorem 1. At the beginning of this 
section, we claimed that Theorem 1 was equivalent to the statement that the subtableau 
has shape approximately described by the region bounded under the graph of the level 
curve {L = k/n 2 } (which in rotated coordinates is given by the curve v = g a {u), where 
a = k/n 2 ). We shall now make precise the sense in which this is true, and see how this 
follows from the fact that g a is the minimizer. 

For a continuous function p : [—y/2/2, y/2/2] — > E, define its supremum norm 

I bl loo = max \p(u)\. 

uG[-V2/2,V2/2] 

Lemma 2. K is continuous in the supremum norm on the space of admissible func- 
tions. 
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1 /-V2/2 /-V2/2 

(9,h) = -- _ / _ s'(<0/i'(i) log |s-i|dsdi 

2 J -y/2/2 J-V2/2 



(10) 



Proof. Consider the symmetric bilinear form 

"V2/2 /-V2/2 
-■v/2/2 J-y/5/2 

defined whenever g and h are almost everywhere differentiable functions on [— y/2/2, y/2/2] 
with bounded derivative. We show that (•, •) is continuous in the supremum norm with 
respect to any of its arguments, when restricted to the set of 1-Lipschitz functions; this 
will imply the lemma, since K(g) = (g,g). Write (|1(J|I more carefully as 

rv / 2/2 

h'(t)log(s-t)dt + I h'(t) log(f - s)dt 



(9,h) 



1 f^ /2 

2 J-V2/2 c \-° 



ds. 



-y/2/2 J s+c 

For s E (-V2/2, y/2/2) which is a point of differentiability of h, integration by parts gives 



h'(t) \og(s-t)dt 



V2/2 



h(t) log(s-t) 



V2/2 



h'(t) \og(t-s)dt = 



t=-V2/2 J-V2/2 t - s 



dt + h(t) log(t - s) 



t=y/2/2 



t=s+e 



y/2/2 



S + € 



hit) 
t-s' 



dt 



+ (h(s - e) - h(s + e))loge 



Ht) 

■%/2/2,s-e]u[s+e,v / 2/2] t — S 



(ft 



e\0 



%h(s), 



where h is the Hilbert transform of h, defined by the principal value integral 



h(s) = - 



h{t) 



dt 



7r Jn.t- a 

(think of h as a function on R which is outside [—y/2/2, y/2/2].) Going back to (fTT)|) . this 
gives 



(9,h) 



2 2 



4* 



-y/2 



V 2 

y/2/2 
■y/2/2 



9 [s)log \ — - s \ ds 
-V2/2 \ 2 / 

V2/2 \ 
5'(s)log \ — + a]da 

-y/2/2 \ 2 / 



g'(s)h(s)ds. 



(11) 
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Now recalling that the Hilbert transform is an isometry on L 2 (M.) (see [531, Theorem 90), 
and using the fact that 



log — ± s ds 
-V2/2 V 2 / 



this implies that for 1-Lipschitz functions g,hi,h 2l 



2 - log 2 



< 1, 



|<5,/ii-/i 2 )| < ||/n — halloo + / _ \hi(s) - h 2 (s)\ds 

z J-V2/2 

< \\h 1 -h 2 \\ OD +2 1 ^^J^ (hi(s)-h 2 (s)y ds^j 

< \\h 1 -h 2 \\ oc + 2 1 /^ ^£{hi(s)-h 2 (s)y ds^j 



1/2 



1/2 



1 - Moo + 2 1/4 J f / ^ (fti(s) - ft 2 (s)) 2 ds 
2 W-V2/2 , 



1/2 



< (l + 2V>I) 



We have another use for Let / be a Lipschitz function on [— v2/2, \/2/2] that 

satisfies /(±\/2/2) = 0. Denote by 



*"[/](*) = / f{t)e- lxt dt 
Jr 

the Fourier transform of a function /. Recall the well-known formulas 

= ». 8 gna:.F[/](a;), 
= i-x-F[/](or), 



Then, by (EJ 



fx(t)f 2 (t)dt = — I F[h]{x)F[f 2 ]{x)dx. 



V2/2 



K(f) = {f,f) = - f(s)f(s)ds 
z J-V 2 / 2 



1 



1 



= I I F[f](x)F[f](x)dx = ~ I \x\ ■ \F[f](x)\ 2 dx. 



(12) 



We note as a lemma an important consequence of this identity which we shall need later 
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Lemma 3. If / is a Lipschitz function with /(±\/2/2) = as above, then K(f) > 0, 
and K (/) = only if / = 0. ■ 



Lemma 3 will be used in the next section to easily deduce uniqueness of the minimizer. 
In fact, Theorem 7 gives all the necessary information to prove a non-quantitative version 
of Theorem 1, i.e. without the rate-of-convergence estimates. However, we can do better, 
by noting that Theorem 7(iii), together with the representation i|12[) . can be used to give 
quantitative estimates for the rate of convergence in Theorem 1. We prove the following 
strengthening of Lemma 3: 

Lemma 4. For every r E (2, 3), there exists a constant c = c(r) > such that for all 
2-Lipschitz functions / : [-a/2/2, a/2/2] R that satisfy f(±y/2/2) = 0, we have 



Proof. Had the power of |a;| in l|12|) been 2, K(f) would have been equal to 1/4 times 
the squared L2-norm of xF[f] (x) = F[f] (x). Having |x| in l|12H invites the conclusion that 



/ of order 1/2. 

To see that this is indeed the case, and to use the full power of such an interpretation 
of K(f), let us recall the corresponding definitions. For a S (0, 1), the fractional derivative 
f( a \x) of order a is defined by 



The integral exists as f(x) is Lipschitz and bounded. Clearly f^ a \x) = for x < — a/2/2. 
Then 



*(/)>c||/n; 



instead we are dealing with the squared L2-norm of /^/^(x), the fractional derivative of 




(13) 




(14) 



where 




|a;| a exp(ia!7r/2), x > 0, 
|x| Q exp(-ia7r/2), x < 0. 



Indeed, setting 



\z\exp(i(l + a)6)), if z = \z\e id , G € (-n,ir), 
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we have 



1 — e~ ixT f l °° 1 — e~ z f°° 1 — e~ T 
— dr — (ix) a / — =- dz = (ix) a / — 7- dr 

T i+ a v ) j o z i+ a v > j q T i+ a 



{ix) a - P T' a e~ T dT = ( ? x)° r(1 a \ 







In particular, for a = 1/2, we get from (|14ll that 

\F[fW](x)\ 2 = \x\-\F[f](x)\ 2 , 
whence, by (|14|l and isometry of the Fourier transform, 

K(f) = \J r \x\- \F[f]{x)\ 2 dx = ^\fW\x)\*. (15) 
The fractional integration operator, inverse to that in i|13[l . is known to be given by 

f{x) = (I a f^)(x), (I a h)(x) := -L- f (x - tr- x h{t)dt. (16) 

As a check, the Fourier transform of the RHS is 

1 f°° 
—Flf^](x) / T a ~ l e- ixT dr = (ix)- a F{f^](x) = F[f](x). 
r (a) Jo 

By Theorem 383 in 15 , for p > 1 and 

1 

< a < -, 



p 1 ~ ap 

I a maps L p into L q , and is bounded. That is, there exists a constant c(p) > such that 

\\I a h\\ q <c( P )\\h\\ p . (17) 

Introduce ip(x) = /^"' ( x )l(-oo V2/2] ( a; )' so tnat ^ ^ s su PP or ted by [— \/2/2, V2/2]. Ac- 
cording to 11611 . 

(/aV)(z) = /(S), X < V2/2. 

So, using Ijl7|) and monotonicity of the L s -averages, we have 

ll/ll, < |M>||, < c(p)M\ p < Cl (p)M\ 2 < c 2 (p)||/W|| 2 , Cl (p) := (V^) 1 ^ 1 /^). 
In light of (|15fl , for a = 1 /2 we obtain then 

~<c 2 (p)K(f), c 2 (p):=- Cl (p) 2 , fp 6(1,2), g =_P_y (18) 
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Let x <E (-V2/2, y/2/2) be such that |/(x )| = ||/||oo- Since / is 2-Lipschitz, 
|/(aOI>||/||oo- 2\x-x \, \x-x Q \<M}hl. 

Then 

/ /•H/lloo/2 \ 2/9 || f ||2(g+l)/9 

ii/ii^( 2 y (ii/iu-2 y )^j = " ( ^ l)2/g , 

so, using l|18[). we conclude that, for an absolute constant c*(p, q) > 0, 

K(f)>c*( P )\\f\\^+^. 

It remains to observe that 

2(g + l) =1 | 2 

can be made arbitrarily close to 2 from above by selecting p sufficiently close to 2 from 
below. This completes the proof. ■ 

Theorem 8. For a Young diagram A whose graph lies within the nxn square, let g\ (u) 
be the rotated coordinate version of the function f\(x) defined in Q. Denote a = k/n 2 . 
Then for all 2 < r < 3, there are constants c = c(r) > 0, C — C(r) > such that for any 
e > and for any n, 



P„ TeT„: max \\g X k - g a ||oo > e < Cexp(3n - ceV). (19) 
V l<fc<n 2 — l T y 

Consequently, with probability subexponentially close to 1, for all k the supnorm distance 
between g A * and g a , (a — k/n 2 ), does not exceed n _1 / 2+l5 , (5 > 0). 

Proof. Let p{m) be the number of partitions of an integer m. It is known that for all 
to, p(m) < exp(7ry / 2m/3) (see 0, Theorem 14.5). Fix n, 1 < k < n 2 — 1, e > 0. Using 
Lemma 1, 

reT n :||<? A , -ffalloo >e^ = P n (r e T n : \ k T = \ 

Ao C □„ of area fc 

||.9Ao - 3a I loo > e 



<p(k) sup exp( -(l + o(l))n 2 (K{g Xo ) + H(a)-log2)). (20) 

Ao C □„ of area A: 
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Let Ao be a diagram contained in □„ of area k, such that \\g\ — <?a||oo > e- Since g\ is 
a-admissible, using Theorem 7 and Lemma 4 we have 

K(g Xo )+H(a)-log2>K(g Xo -g a )>c(r)\\ 9Xo -g a \\ r 00 >c(r)e r . 

Combining this with 12U|) and with the above remark on the number of partitions of an 
integer gives that for n larger than some absolute initial bound, 



TeT n : || 5a * -g a \\oo > ej < exp(2.8\/an - cn 2 e r ). 
Taking the union bound over all 1 < k < n 2 — 1 gives <(TT)|) . ■ 

Lemma 5. For each (x, y) £ (0, 1) x (0, 1), let (u, v) be their rotated coordinates as in 
(Q. Let ao — L(x,y), so that |u| < y / 2oo(l — ao) and v — g aa {u). There exist absolute 
constants C\ , c 2 > such that if we set 

a(x,y) = min(.Ty, (1 - x)(l - y)), 

d(x,y) = cx-y/a(x,y), A(x, y) = c 2 cr 2 (a;, y), 

we will have that for all < a < 1 and <5 < A(ie, y), if — .g QQ (u)| < 5 ■ d(x, y) then 

\a — ao\ < S. 

Proof. Since g a (u) increases with a, it suffices to prove existence of two absolute con- 
stants 71 , 72 > such that 

\g a {u) - g aa (u)\ > ~fi<T 1/2 (x,y)\a -a \, if \a>- a \ < j 2 cr(x,y). 

Because of the symmetry property g\- a {u) = — g a (u), we may assume that x + y < 1, 
or equivalcntly that ao < 1/2. 

To prove the above claim, we note the following inequalities. Notice first that 



y/2a (l - a Q ) > v =4> a > 



1 - VT^2 



v- 



2 

Likewise, a^ that corresponds to the lowest point (u,u) is given by 



a y ' — 



and we see that 



~0 2 2 2 



i \ vl — 2u 2 — vl — 2v° v~ — w v~ — u" 

a - a { ~> > = ; ; > = xy. (21) 

2 VI - 2w 2 + VI - 2v 2 " 2 1 ; 
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(|21[1 says that decreasing ao by 2o2/o gives us a feasible a, for which (u,g a (u)) lies between 
(u, v) and the lowest point (u, u), such that u < y2o(l — a). 

Let us estimate from above g a (u) for a <E [cr ao]- From l|t)()|) it follows that 



dg a (u)/da 



> c 



for some absolute constant c > 0. (Indeed, 2a(l — a) = (3 2 (a) > g a (u) 2 .) Integrating from 
a G [or },a!o] and exponentiating, we obtain 



g ao (u) + y/g ao (u) 2 -u 2 
g a (u) + ^fg a (u) 2 - u 2 



or equivalently 



Consequently 



g a (u) - \/g a {u) 2 - u 2 



l ao {u) - \Jg ao {u) 2 - u 2 



> exp(c(a - ")), 



> exp(c(ao — a)). 



\/g a {u) 2 - g ao ( u ) 2 < cosh(c(a - a))\/g aQ {u) 2 - u 2 - sinh(c(a - a))g ao (u), 



or 



g a (u) 2 < 



cosh(c(a - a))g ao (u) - sinh(c(a - a ))^g ao (u) 2 - v? 



so that 



Ja(u) < cosh(c(a - a))g ao (u) - sinh(c(a - a ))V 9a () (u) 2 - u 2 



Consequently, for some constants a > 0, 

g a (u)-g ao (u) < -c 3 (a - a)[(v 2 - u 2 ) 1/2 - c 4 (a - a)v] 
= -C5 (a - a) \{xy) 1/2 - c 6 (a - a) (x + y)] 
< -c 7 (a - a)(xy) 1/2 , 



provided that 



a - a < c 8 



(xy) 



1/2 



x + y 



From (EB we know that we can go below ao by xy at least. Pick p = min(l, Cs/3); then 
the last inequality holds for — a < pxy, and we have 



g a (u) - g ao (u) < -c 7 (a - a){xy) 1/2 , a € [a - pxy,a ] 

Now for «o < ol < 1/2 we know that 

dg a (u) 



(22) 



da 



> co, \/v 2 - u 2 = c w {xy) 1/2 , 
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so that 

g a {u) - g ao (u) > c w (xy) 1/2 (a - ao). 
By symmetry, for l/2<a<l — ao, 

9a(u) - 5i-ao < -cio((l - «o) - «)((! - sc)(l - 2/)) 1/2 , 
and, for 1 — ao < < 1 ~ ceo + p(l — — y)i 

~g a {u) - g Q0 (u) > c 7 (a - (1 - tt0 ))((l - z)(l - y)) 1 ' 2 . (25) 
The inequalities J23), prove the claim with 71 = min{c7, cio} and 72 = p. 



(23) 
(24) 



Proof of Theorem l(ii). We now prove Theorem 1(h), the part of Theorem 1 
that deals with the interior of the square. The treatment of the boundary of the square is 
more delicate and is deferred to section 4, being essentially equivalent to Theorem 3. 
Fix 1 < i,j < n such that 

mm(ij,(n-i)(n- j)) >n 3/2+e . (26) 

Let (u,v) be the rotated coordinates corresponding to (x, y) = (i/n, j /n). Let ao = 
L(i/n, j/n), so that v = g ao (u) and |w| < -\/2ao(l — ao). For each tableau T — (iij)i<i,j< n G 
T n let kx = tij, and let j3r = kx/n 2 . Note that kx is an integer representing the smallest 
s such that X T contains the box (i,j). This implies that 

, \ V2 
v < g x k T (u) < v H , 



V2 

j x k T (u) - g ao (u) < . (27) 

Apply Lemma 5 with (x,y) = (i/n, j/n) and S = n~^~ e '^ 2 . Note that because of 
for n large we have S < A(x, y) as required. Then, making use of Theorem 8 we get 



TeT n 



1 



2 



> 



1 



n \n n 

(by Lemma 5) < 



< 



reT„: |#r-ao| > 



n (l-e)/2 J ~ ~„^^„.^ -u, - (i_ e)/2 
P„(T G T„ : ^(u) - g ao (u)\ > ^fcff ) 



(by C3 

(for n large enough, by l|26|l) < 



5 fc T (u) -g(3 T (u) 



(j 1 e % : g^ (u) - gp T (u) 



> 



> 



d(i/n,j/n) _ V2\ 
n(i-0/2 n J 

d(i/n,j/nY 
2n( 1 - £ )/ 2 



(by Theorem 8 with r = 2 + e) < C exp 3n — cn 



2 / d(i/n,j/n) 



2n( 1 - £ )/ 2 



2+e N 



(for n large, by < C"exp(-c' 
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Taking the union bound over all 1 < i,j < n satisfying l|26|) gives the result. 



3 Solution of the variational problem 
3.1 Preliminaries 

In this section, we prove Theorem 7. We actually derive the explicit formula for the 
minimizer using methods of the calculus of variations and the theory of singular (Cauchy- 
type) integral equations. Our derivation makes only one a priori assumption (obtained 
by educated guesswork and later verified by computer simulations) on the graphical form 
that the minimizer would take, and so is in a sense more systematic than the analogous 
treatments in the fundamental papers |20|. |35|. |36| . where the solutions are brilliantly 
guessed using the properties of the Hilbert transform. We believe that our technique may 
prove useful in the treatment of similar problems in the future. 

First, observe that because of symmetry, we need only treat the case a < 1/2; the 
mapping g — > y/2 — g takes the set of a-admissible functions bijectively onto the set of 
(1 — a)-admissible functions, and has the property that K{-\/2 — g) — K(g). 

Next, observe that for a = 1/2 the assertion is immediate, because of Lemma 3. 

We prove another fact that follows from general considerations, before turning to the 
derivation of the minimizer. 



Lemma 6. For any < a < 1, the functional K has a unique a-admissible minimizer. 

Proof. The functional K is continuous on the space of a-admissible functions, and is 
bounded below by Lemma 3. By the Arzela-Ascoli theorem, the space of a-admissible 
functions is compact in the topology induced by the supremum norm (since the admissible 
functions are uniformly bounded and equicontinuous). Therefore K has a minimizer. To 
prove that the minimizer is unique, let h\ and hi be two distinct a-admissible minimizers. 
Then h — (/ii+/ia)/2 is also an a-admissible function, and g = (hi— h^jl ^0, g(±-\/2/2) = 
0. So, using the parallelogram identity and Lemma 3, 

K(h) = \K[h{) + \K{h 2 ) - K{g) < . min K(h), 

Z Z n. is a-admissiblc 

a contradiction. ■ 
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3.2 The derivation 



We now proceed with the derivation of the minimizer, which we shall denote h = h a . 
The dependence on a will be suppressed except where it is required. For the rest of this 
section, a will be a fixed value in (0, 1/2), unless stated otherwise. 

First, note that, under the condition h(±-\/2/2) = \/2/2, the a-condition J_^/ 2 Q l ( u ^ 
\u\)du = a is equivalent to 



V2/2 1 

uh'(u)du = a . (28) 

V2/2 2 

We now formulate a sufficient condition for h to be a minimizer. It is based on a stan- 
dard recipe of the calculus of variations, the Lagrange formalism. We form the Lagrange 
function 

rV2/2 

£(h, A) = K(h) — A / _ uti(u)du 
J-V2/2 

and require that, for some A, h a be a local minimum point of C(h, A) in the convex set 
of functions h subject to all the restrictions except the a-condition J2HJ- To be sure, we 
ought to include into the function a term A' times the integral of h! , since h must meet 
another constraint 

\/2/2 

h'(u)du = 0. (29) 

V2/2 

We chose not to, since - in the square case - even without this constraint h'(u) will turn 
out to be odd anyway. Since C(h, A) depends explicitly on h! alone, we get the equations 
for the sufficient condition in a simple-minded manner, by taking partial derivatives of C 
with respect to h'(s), s G (— v2/2, V2/2), and paying attention only to the constraint 
— 1 < h'(s) < 1. The resulting "complementary slackness" conditions are 

r V2/2 f =0, if - K h'(s) < 1, 

w(s):=- _ h'(t)log\s-t\dt- As is<^ > 0, ifh'(s) = -l, (30) 



<0, iih'(s) = l. 



Lemma 7. If h is an a-admissible function that, for some Agl, satisfies (|3U|) for all 

s e (-V2/2, y/2/2) for which h'(s) is defined, then h is a minimizer. 

Proof. If g is a 1-Lipschitz function on [s/2/2, V2/2], then (j20| implies that (g'(s) - 
h'(s))w(s) > for all s for which this is defined, so 

y/2/2 /-V2/2 

g (s)w(s)ds > / h! (s)w(s)ds. 



-V2/2 J-V2/2 
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If g is a-admissible, by l|28[) this can be written as 



1 



2(h,g)+a-~ = 2{h,g)-X sg'(s)ds 



V2/2 



2 J-V2/2 

> 2(h,h)-X _ sh'(s)ds = 2(h,h) +a- -, 
J-V2/2 2 

which shows that 

(M > (M>- 

Therefore, by Lemma 3 applied to the function g — h, 

(g, g) = (h, h) + 2{h,g-h) + {g-h,g-h)> (h, h) , 
so h is a minimizer. ■ 

We are about to prove part (i) of Theorem 7, namely that h — g a is the minimizer. 
Assuming this, note that in the above proof we actually showed that 

(9,9) > (9a,9a) + (g~9a,g- 9a), 

which is precisely the claim of part (iii) of Theorem 7. So it remains to prove parts (i) 
and (ii). 

Our challenge now is to determine an admissible h that meets the conditions . 
Now look at Figure 1(c) with your head tilted 45 degrees to the right. Based on the shape 
of the level curves, we make the following assumption: For some (3 = (3(a) & (0, \/2/2), 

= -1, if - V2/2 < s < -13, 
h'(s) is <j e(-l,l), i£-p<s<P, (31) 
= +1, if (3 < s < y/2/2. 

Substituting this into (|5U|l gives that for —f3 < s < (3, 

rP rP rV2/2 

- h'(t)log\s-t\dt = Xs- _ \og{s-t)dt+ \og(t-s)dt 

J-P Jp 

= Xs + (V2/2 - s) log(V2/2 - s) - (V2/2 + s) log(V2/2 + s) 

+ (f3 + s) log(/3 + e) - {ft - s) log(/3 - s) (32) 

Assume that h'(s) is continuously differentiable on (—(3,(3). Differentiate l)32|) . to obtain 

13 h'(t) B 2 - s 2 

IStldt = X + log^-, (33) 
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where the left-hand side is a principal value integral. 

In the theory of integral equations this is known as an airfoil equation. Solving it is 
tantamount to inverting a Hilbert transform on a finite interval. Fortunately for us, it 
can be solved! The following theorem appears in ^U], Section 3.2, p. 74. (See also |28|. 
Section 9.5.2.) 

Theorem 9. The general solution of the airfoil equation 

i f 1 g(y) 



■dx = f(x), \x\ < 1, 

■n J-iV-x 

with the integral understood in the principal value sense, and fix) satisfying a Holder 
condition, is given by 



*(*) = * f Vrr¥m dy + -^=. (34) 



Applying Theorem 9 to l|33[) . we get the equation 

*-<■> = Twhw ll {02 - 12,1,2 ( A + log t£) + 



(/J2 _ s 2)l/2 



(35) 



2 

Here the integral is again in the sense of principal value, and the equation must hold for 
some value of c. 

We evaluate the integral in l|35|l . Consider the contribution of the A-term first. Sub- 
stituting t — [3 sin a; and later u = tanx/2, we get 

dt = (3 I — : — dx 



s — t J-tt/2 s /P ~ sina; 

li 

(s//3 + sinx 

/2 

2/3 (l - (s//3) 2 ) f 1 du 



r /2 r /2 dr 

(3 (s / (3 + sin x)dx + P (I ~(s/f3) 2 ) — — 

J-tt/2 J-k/2 s/p-smx 



s/(3 J-t u 2 - 2(/3/s)u+ 1' 

For |s| < /3, the denominator in the last integral has two real roots, u\ G (—1, 1) and 
U2 ^ (— 1,1)- A simple computation shows that the principal value of this integral at 
u = u\ is zero. So 

r (/3 2_ £2)1/2 



-a s-t 



-dt = irs, s€(-f3,f3). (36) 



Turn to the log-part of the integral in i|3"5|) . Substituting t — t/3, s = V\(3, (2/3 2 ) 1 = v 2 , 
we see that 



J-3 



(/3 2 _ £2)1/2 01 _ t 2 



lo S ^ — T^dt = /3[I(s//3, 72/(2/3)) - I(-s/p, V2/(2/3))], (37) 



s-t 
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where 



j(e.7) 



1 q~v) lf \ 1 + 17, 

— ; log — — dr;, 

-l 4 — ^7 7 + ?7 



ee[-l,i],7>i, 



is evaluated in the following lemma. 



Lemma 8. 



J(£,7)=7T 



1 - 7 + V7 2 - 1 - £ log (7 + vV " l) 



2 a/!- tan 



-1 (7 

V(7 + l)(l + 



(38) 



Proof. Notice that /(£, 1) = 0, and, for x > 1, 
0I(£,a;) _ /- 1 (1 - 77 2 ) 1 / 2 



dx 



1 — »7) (a: + »7) 



dij 



1 (1-V)V2 ,1 (1-^)1/2 

7 «?? + / an 

•i ?~»7 



1 x + 77 



< 1 ^(1-^)1/2 



x + £ + x + T] 



-drj, 



(39) 



see ||5d]) . Substituting 77 = sini, (i £ [— 7t/2, 7r/2]), and then t = 2tan 1 u, (u £ [—1,1]), 
we evaluate 



1 (1 _ ^2)1/2 

■1 a; + ?7 



dry = (xt + cost)^ + (1 - x 2 ) 



^ 2 dt 



x + sint 



ttx + 2(1 - x 2 ) 



du 



! x(l + w 2 ) + 2it 



= ttx - 2(x 2 - l) 1 / 2 

= 7r(x — (.T 2 - l) 1 / 2 ) 

Combining this with i|39|) . we obtain 

9a; 



tan 



a; + 1 
x- 1 



tan 



1 >~1 
x + 1 



(40) 



^(x 2 - l) 1 / 2 

x~+£ 



We integrate this equation from x = ltox = 7>l, and use the substitutions x — coshi, 
t £ [0,t ], with 

to — arccosh7 = log ^7 + (7 2 — 1) 
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and then u = e*, u £ [l,uo], with 

u = e to =7+(7 2 -l) 1/2 - 

We have 



=-7r(7-l)+7r 

= — 7r(7 — 1) + 7T 

Here 



sinh 2 t 
cosh i + £ 



(sinht-^)|*°+ 2(^-1) 



o w 2 + 2£u + 1 



- 1 



(sinht - £t)|*° = ( 7 2 - 1)V2 _ ^ log ( 7 + ( 7 2 
and the last integral equals 

, K-i)(i-£ 2 ) 1/2 



: tan 



(1_ £2)1/2 
1 

: tan 



: tan 



uo-i i + e 



l-^ + (uo +0(1 + 

i 



VT^ 2 uo + 1 V 1 - £ v/W 2 
Combining gTJ, g2J>, 63 gives 

Now from © and gSJl we get 

r- c I 

h'(s) 



( /3 2_ s 2)l/2 n (02 _ s 2)l/2 



A -2 log 



\/2/3 



H — tan 

7T \ 



-i /(7-l)(l + 



tan 



. (7 + l)(l"0 
with £ = s/p, 7 = 72/(2/3), or, after some simplification, 



h'(s) 



A - 2 log 



/(7- 


1)(1- 




/(74 


■i)(n 




14 


Vi- 


2/? 2 


og — 







(41) 
(42) 



V(7+l)(l + 1 ' 



(44) 



(/3 2 -s 2 )V2 tt(/3 s - s 2 )V2 
I 2 tan- 1 ( 1 - 2 ^) 1/2s 

+ ^ tan ( /3 2_ s 2 ) l/2 ■ 

We now observe that the only values of c and A for which the right-hand side is bounded 
as s /* P, s \ — P, and therefore has a chance of being the derivative of an a-admissible 
function, are 



c = 0, A = 2 log 



i + yr^w 

V2P 



(45) 
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Therefore we get 

»'<•> = <«) 

Note that h'(s) e (—1, 1). We have determined h'(s), except the value of j3 = (3{a) such 
that h is a-admissible, i.e., satisfies (123 ■ Rewrite Egl as 



/3 

sh'(s)ds = a-{3 2 . (47) 

-/3 

Besides evaluating this last integral, to compute h(s) explicitly we will need f^_g h'(u)du. 
To this end, integrating the first arctangent-of-radical function in (|44|l on the interval 
[-U], (£e (-1,1]), we get 



l 



( 7 + l)(l-f?) 



V(7+l)(l-r?) 2 y.j^-^^/rr^a v ' 

Substituting in the last integral 77 = sint, and then u = tani, we transform it into 

- to - 7; + 7 / r- to = sin £ 

2 7-^/2 7 -suit 



7T 



mi 



- t0 -2 +2 L l + u*-2u h K = tant /2] 

t g I 27 M °' 7 ' 1 I tan- 1 1+7_1 



2 Vr^V v/i^t^ '" V^T^ 



n 27 _ 1 /l + u /7— 1 



-i -- + ^^tan- 1 ^ ^ ; (49) 



here 



1 + m _ 1 + tanio/2 _ 1 + sint _ 1 + £ _ + ^ 
1 - w ~ 1 - tan t Q /2 ~ cos t ~ ^/l - £ 2 ~~ y 1 - T 

From gSJl, © we obtain 



? tan- J^P^rd, 
-1 V(7 + l)(l-^) 



(50) 



» (*"' £+ 5)- (51) 
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In a similar fashion 



,tan- f^ + _f^ r, = J(l - 7 2 + 7V^D, (52) 

and the integral in the negative arctangent in 144f> is obviously given by (|52[1 as well. Using 
(|44|l and (|52() . we see that the a-condition (|47|l is equivalent to 

/? 2 ( 7 2 - 7v / 7 2 - 1) = a 1 - 2a = a/1 " 2/3 2 , 

the latter being possible only if a < 1/2. In that case 



= v/2a(l - a). (53) 

Consequently, see 



A = l 0g L^ 7 (54) 
a 

and, see (|4*^|) . 

, s € (- V2a(l - a), a/MI - a)). (55) 
Furthermore, denoting the integral in (|51J) by J(£, 7), we easily get 



/i'(s) = - tan 1 v — 

71- y^/2a(l-a)- 



h(s) = + h'(t)dt = 0(1 + J(£, 7) + J(-£, 7 ) - J(l, 7)) 

J-/3 



= 2 -s tan^ - (1 = 2 ^ + ^ tan^ Eg Eg , (56 ) 

7T ^ ^2a(l - a) - s 2 y 7T ^ l-2a J 

We have derived a formula for a candidate minimizer, which we now recognize as the 
function g a that we defined in section 2. To be sure, this function was determined so as 
to meet the ramifications of some of the constraints. However, looking at 1)55(1 . we see 
that —1 < h'(s) < 1 for s £ (—0,0), so h is indeed 1-Lipschitz, even though so far we 
haven't paid attention to this constraint! Furthermore, since h'(s) is odd, the constraint 
(|29|l is met automatically, and it is the reason why we were able to satisfy the boundary 
constraints h(—y/2/2) = h(s/2/2) — v / 2/2. Also, we determined from the requirement 
that h should satisfy l|47|l . which under these boundary conditions is equivalent to the 
a-condition. We conclude that, at the very least, g a meets all the constraints, thus is 
a-admissible. 

By Lemma 7, to prove that g a is the minimizer, it only remains to prove that g a 
satisfies the conditions (J3UJ). By (SHI, w'(s) = for \s\ < 0. And w(0) — as h'(t) is odd. 
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So w(s) = for \s\ < (3, hence the first condition in (|30fl is met. As for the remaining 
conditions, by (anti)symmetry, it suffices to check, say, the third condition, namely that 



/•a/2/2 

F(s,a):=- _ g' a (t)log\s-t\dt- X(a)s < 0, /3(a) < s < a/2/2. 

^-\/2/2 



Fix < s < a/2/2, and let a = (1 - a/1 - 2s 2 )/2, so that /3(a) = s. Clearly, because 
of the first condition in H30|) . F(s,&) = 0. To finish the proof, we will now show that 
dF(s,a)/da > for < a < a. By O, 

9F M- ^^^^logl,-^-^. (57) 



dp df3 01 1 9/3 

Using lfi?)|l and simplifying gives 

dg'Jt) = 2 t 

9/3 tt/3(1 - 2/3 2 )V2 Qga_t2)i/2- 

Since /3'(a) = (1 - 2/3 2 ) 1 / 2 //3, J57J) becomes 

dF(s,a) 2 / >/3 £log|s-i| 



9a tt/3 2 (/3 2 -t 2 )!/2 T (i_a) a ' 
Here the integral equals 



09a-*3)i/2_i og | fl _ t | 



-/3 J-P s ~ l 



(p2 _ +2\l/2 



see l|5B). Therefore 



_ _ilf a _^ 2 « 2 W 2 



( a -( s a-/3 a )Va) + 



9a /3 2 V / (1 — a)a 

2 ^ +-^( S 2 -/3 2 ) 1 / 2 



(1 - a)a /3 2 J (3 

3.3 Direct computation of K(g a ) 

Our next goal in this section is to show that K{g a ) = —H(a) + log 2. There are two 
ways to do this. First, looking at the proof of Theorem 8, we see that we may repeat the 
arguments of that proof (without assuming the value of K(g a ) as in that proof) to deduce 
that the value M a of K(g a ) + H{a) — log 2 must be 0. For, if it were greater than 0, then, 
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denoting k = [an 2 \ , we would have 

i=P„(TeT n ) = J2 

Ao of area k 

< p(n 2 ) exp (-(1 + o{l))n 2 M k/n 2) ► 



(since M a is obviously continuous in a.) On the other hand, if M a < 0, then for some 
sufficiently large n, we would have for some diagram Ao of area [cm 2 J contained in □„, 
that K(g\ ) + H(a) — log 2 < (take a diagram for which g\ approximates g a , and use 
Lemma 2). But this again implies a contradiction: 

1 > P„ (T e T n : AL? n2j = Ao) = exp (-(1 + o{l))n 2 (K(g Xo ) + H{a) - log 2)) > 1. 

These last remarks notwithstanding, we find it worthwhile to compute K(g a ) directly, 
if only to thoroughly test our derivation of g a , and to show that all the integrals involved 
can be evaluated explicitly. 

For h = g a , rewrite l|2*S|) as 



%/2/2 

-V2/2 



u(h (u) — sgnu)du = a. 



Using this, multiply both sides of PU)> by (h'(s) — sgns) and integrate, obtaining 

,72/2 



K(h) 



Xa 1 
~2 ~ 2 



-V2/2 



h'{t) 



2flog|t| - (* + v / 2/2)log|t + V2/2| 



dt, 



-(t-V2/2)log|t-V2/2| 

where we found before that A = log((l — a) /a). Denote 

S(t) = 2t log \t\ - (t + V2/2) log \t + V2/2| - (t - V2/2) log \t - \/2/2|, 

and set 



(58) 



V2/2 



-V2/2 



h'(t)S(t)dt, 



so that K(h) = -Xa/2 - K t (h)/2. Just like JS7JI, 



^l(feq) 

a/? 



\/2/2 



72/2 9/3 
•/3 



S(t)tft 



tt/3(1-2/3 2 ) 1 /2 J 



(/3 2 - i 2 ) 1 / 2 



7T/3(1 - 2/3 2 )V2 ,/_^ (/3 2 _ (2)1/3 

2 log |t| - log |t + V2/2| - log |t - V2/2| 



dt. (59) 
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Denote 



so that 



E(s,/3) = / (P 2 - t 2 ) 1 / 2 log \t-s\dt, 



dK g { p a) = 2E(fl, /?) - E(-V2/2, (3) - E(V2/2, (3). 



(60) 



By and 

dE(s,P) _ f (/3 2 -t 2 ) 1/2 



ds 



Therefore 



s - t 



dt 





TTS \s\ < (3, 

7r(sgn S ) - (s 2 - /? 2 ) 1/2 ) , (3<\s\< V2/2. 



(61) 



2E(0,/3) = E(f3,/3)+E(-/3,j3)+ir / sds + ir j sds = E{f3, 13) + E(-fi, (3) - -Kf3 z . (62) 



Likewise 



E(-V2/2,0)+E(y/2/2,0) = E(-0,0)+E(p,0)+2n / (s - (s 2 - /3 2 ) 1/2 ) ds, (63) 



V2/2 



where 

«. a/2/2 
/3 



( s 2 _ ^2)1/2^ = 1 J s(s 2 _ ^2)1/2 _ 01 lQg ^ + (s 2 _ ^2)1/2^ 

1 (\-2a 



2 V 2 



a(l — a) loj 



V2/2 
1 - a 



(64) 



So, using /? = (2a(l - a)) 1/2 , 



E(-V2/2, /3) + E(V2/2, (3) = E(-0, 13) + E(J3, f3)+ir l^-f3 2 + a + a(l - a) log 
and, combining this relation with l|62|) . we simpify (|60|l to 



1 — a 



So, by P|) 



— — — = — 7r a + a(l - a) log 

op \ a 



dKxiK) 8K x {h a ) (l-2(3 2 ) 1 / 2 1 1-q 

log . 



da 



8(3 







1 - a 
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Since h' a = at a — 1/2, we have Ki{h) = at a = 1/2. Hence 
Ki{K)= I (— ^+log^— ^<te = -log(l-a)-21og2 

Jl/2\1-X X ) 

— (1 — a) log(l — a) — a log a, (65) 

which gives finally for K{h a ) 

K(h) = a log a + (1 - a) log(l - a) + log 2 = -H(a) + log 2. 
The proof of Theorem 7 is complete. ■ 



3.4 The parametric family g a 

The minimality proof in section 3.2 relied on the possibility to consider simultaneously 
the whole family of variational problems, and thus to differentiate the minimizcr g a with 
respect to a. Moreover, to reveal a little secret, we anticipated the formula (|54(l for 
the Lagrange multiplier A. According to a general (semiformal) recipe of the calculus of 
variations (more specifically, mathematical programming), we knew that this A, dual to 
the a-condition, should be equal to dK(g a )/da, which we have proved to be correct. The 
advantages of this approach of varying the parameter a go even deeper than that. It will 
turn out that the partial derivative of the minimizer g a { ) with respect to a is the key to 
the distributional properties of the random tableau. Using the formula for the minimizer, 
we compute easily that 



dgjMj _ \ -Jlail-a) < |u| < 

For each a, direct integration reveals that dg a (u) / da is a probability density function, i.e. 

dg a {u) 



du = l. 

V2/2 oa 

(In fact, it is the density of the semicircle distribution, and it will play a prominent role 
later. See sections 5, 8.1, 8.2.) This observation is in perfect harmony with the fact that 
g a satisfies the a-condition, thus providing a partial check of our computations. Indeed 

V2/2 r\/2/2 

{g a {u) - \u\) du = I (g a (u) - g (u)) du 

V2/2 J-V2/2 



-V2/2 \./o ds J J Q \J-^2/2 ds 
1 ds = a. 
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Had we been presented with the minimizer g a "out of the blue" , this would have been the 
least computational way to prove its a-admissibility. 



4 The boundary of the square 
4.1 Proof of Theorem 3 

In this section, we prove Theorem 3. As was remarked in section 1.3, the RSK corre- 
spondence induces a correspondence between minimal Erdos-Szekeres permutations ir of 
1, 2, . . . , n 2 and pairs T\, T 2 € T n of square tableaux. By the well known result of Schen- 
sted in this correspondence the length of the longest increasing subsequence in 
7r(l), 7r(2), . . . , 7r(fc) is equal to the length Xj, (1) of the first row of Xj, . So the distribution 
of l n ^k under a uniform random choice of minimal Erdos-Szekeres permutation 7r is equal 
to the distribution of the length of the first row of Aj> in a uniform random square tableau 
T € %i- Denoting for the remainder of this section a = a(k) = k/n 2 , we can therefore 
reformulate Theorem 3 as stating that 



max P n T G T n 

at <k/n 2 <l/2 V 



Xi(l) - 2Ja{l - a)n > aJ 2 u)(n)n) ► 0. (67) 



Theorem 8 looks as if it might imply l|67|l . In fact, it only implies a lower bound on A^(l). 
The reason is that can be very close in the supremum norm to g a (as is known to 
happen with high probability by Theorem 8), while n~ 1 A^(l) might still be much larger 
than 2i/a(l - a) (see (jSHl below). 

Lemma 9. Let a = n,- 2 / 3+e , 8 = n^ 1 / 3 ^-^, e S (0, 2/3). Then 

P n I T € T n : min (A£(l) - 2 v / a(l - a)n) < -Sn = 0(^ b ) (68) 

\ ao<a<l/2 / 

for every b > 0. 

Proof. We use the notation of Theorem 8. The length of the first row A^(l) can be 
extracted from the rotated coordinate graph g x t^ using the following relation: 

-A£(l) = V2 inf ju € [0, V2/2] : g X k_(u) = u} . (69) 



It follows from l|55l) that, uniformly for a € [a , 1/2] and |m| < ^2a(l — a), 



2 + , VMi-") 



\dg a {u)/du-l\ = - tan" 1 V v - > c( ^2a(l - a) - |tt|), 

7r (1 — 2a)|w| 
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c > being an absolute constant. Consequently, for a € [ao, 1/2], 

g a (V2a(l - a) - S) - (v/2a(l - a) - S) > c5 3 ' 2 . 
So if T £ T n has the property that, for some k in question, 



A|(l) - 2\/a(l - a)n < -<Sn, 

then by ffity. 

ll5A^-5a||oo > sup{|.g A fc (u) -3 a (u)| : \/2a(l - a) - 5 < m < y^iap - a)} 



J n ( T £ T n : mm 

ao<c 



= sup{g Q (u) - u : y/2a{l - a) - 8 < u < y/2a{\ - a)} > cS 3/2 . 
So, by Theorem 8 with e := c<5 3 / 2 , 

min (A^(l) - 2y/a{l - a)n) < -Sn 

Ca<l/2 J 

?a x , I Iff** - 9a | |oo > c<5 3/2 ) 

< exp(3n - cn 2 ,5 3r/2 ) < exp(3n - cn 2- ^ 1- ^ 2 ) ► 0, 

n — >oo 

provided that we choose a feasible r, i. e. r € (2, 3), such that r < 2(1 — e) _1 . ■ 

To prove the upper bound and thus conclude the proof of Theorem 3, it suffices to 
prove an upper bound for the expected value of A*, namely that, for ao < a < 1/2, 



<P n [TeT„: max 

Qo<a< ] 



[A£(l)] < 2Va(l-a)n + 0(aJ /2 7i), 



(70) 



where E„ denotes expectation with respect to the probability measure P n . Indeed, choosing 
to(n) — > oo however slowly, we bound 

P„ (T G T n : A^(l) > 2^(1 - + al /2 uj(n)nj 

(by Markov's inequality) < (aj|/ 2 u;(n)n) _:1 'E n max(0, Ay-(l) — 2yJ a{\ — a)n) 
(by Lemma 9, for any b > 0) < (aJ^nM -1 (e„ a£(1) - 2^(1 - a) + <5n + 0(n 1 " fc )) 



O((a /2 n + ^/(aJ^wHn)) = 0{uj{n)- 1 ). 



1/2 



Write 



A* (!)=£> 



J; 



where = Aj.(l) — A^ (1) = indicator of the event that A^ is obtained from A^ 1 by 
adding a box to the first row. Let p n j = E„(7 n j). 
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Lemma 10. In the notation of Lemma 9, as n — ► oo, 

Pnj < " 2 ~ 2j + 0(Sn(n 2 2j + l)" 1 ), 
n^j{n z - j) 

uniformly for ao < j/n 2 < 1/2. 

Proof. Let y n> j be the set of Young diagrams of area j contained in the n x n square. 
For a diagram A € 3^nj , denote by next (A) the diagram obtained from A by adding a box 
to the first row. Then, conditioning I n j on the shape A jT , we write 

A^ = nextCA^ 1 )) = £ rf(A)rf(D„ \ next(A)) 

^ d(next(A))d(D„ \ next (A)) d(X) 

TO d(next(A)) 

This is nearly an average over 3^i,j with respect to the measure {7J|; in fact, slightly less, 
since not any A' 6 34.,j is of the form next(A) for some A € It follows from the 

convexity of the function x —> x 2 that 

2 d(next(A))d(D„ \ next (A)) / d(X) 



d(D n ) \d(next(X)) 

d(X)d(D n \ A) d(X)d(a n \ next(A)) 

d(D n ) ' d(next(A))d(D„ \ A) " 1 j 

AfcJ'n.J— 1 



We now note the amusing identity 

d(A)d(D„ \ next(A)) n 2 - A(l) : 



d(next(A))d(D„ \ A) j^-j + l)' 



(A G 1) (72) 



which follows from writing out the hook products for d(-) in J^Jl and observing cancellation 
of almost all the factors - see Figure 5. Here is a proof of (|72|l . Clearly the only hook 
lengths influenced by this operation are of the cells in the first row and the (A(l) + l)-th 
column. In particular, 

d(X) = 1 ^ X(l)-i + l + X'{i) _ 
d(next(A)) ~ j ^ A(l) - % + A'(i) ' 

here X'(i) is the number of cells in the i-th. column of A. Clearly the fraction factors 
"telescope" on each subinterval of [1,A(1)] where X'(i) is constant. Let [ii,^] be such 
a (maximal) subinterval. Maximality implies that (^A'^)) is a corner of A, and that 
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(A'(ii) + 1, i\) is a corner of □„ \ next(A). Then 



Ml) 

n 



A(l) + l + X'(i) _ A(l) - ii + 1 + A'(ti) _ ^□„\„ cx t(A)(A'( i i) + 1, A(l) + 1) 



A(l)-» + A'(i) 



X(l)-i 2 + X'(i 2 ) 



h\(l,u 2 ) 



where, say, h\(u, v) denotes the hook length for a cell (u,v) S A. Multiplying these 
fractions for all such subintervals i 2 ], we get 

-l 

d(A) 1 



d(next(A)) j 



n /(«>«) • n 9(u,v) 

v (ii,i;)Gcorners(A) / \ (u,-u)Gcorners(D T i.\A) 



(73) 



Here corners(/z) is the corner set of a diagram /i; f(u,v) is the hook length of a cell in 
the first row of A whose vertical leg ends at the corner (it, v) € corners(A); g(u,v) is the 
hook length of a cell in □„ \ next (A) from the (A(l) + l)-th column whose horizontal arm 
ends at the corner (it, v) £ corners(D„ \ A). Next, considering separately the first row cells 
(l,w), v > A(l), the top A'(l) cells in the (A(l) + l)-th column, and finally the bottom 
n — A'(l) cells in that column, we obtain 



d(D„\next(A)) n-A(l) 



A'(l) 



i n 



A(l)-A(fc) 



d(n n \ A) 

Here, analogously to the d(A)/d(next(A)) case 
Y(i) 



A(l) + n 

"' ~ 1 1 k=2 A ^ ~ A ( fc ) ' k 1 A| 1 : + A '(!) ' 



(74) 



1 



n 



X(l)-X(k) + k 



A(l) + A'(l) 11 A(l) - A(fc) + fc - 1 



. (it,t;)£corners(A) 

Putting CZ3J, O, CSJ together gives (|72)> . 
Combining H71|l and (|72|l gives that 



< 7 < E„ 




n 9{u,v) 

(ii,u)Gcorners(D ri \A) 



(75) 



-4 _1 (i) 2 



(76) 



By Lemma 9, we may write 



(8 = n-^ 1 -')/ 3 ), for all j/n 2 e [a , 1/2]. So, using E^[A^T 1 (1)] < E[(A^T 1 (1)) 21 



p 2 < (» a -2j) 



4<5 



••/' : " J i) \ .7) 
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Figure 5: Illustration of ((72} for A = (6,6,6,6,5,5,5,3,3,2): The numbers 
in the cells are the hook lengths before and after the new cell is added. 



or, using (1 + z) 1 / 2 < 1 + z/2 for j < n 2 /2, 



Pn,j < 



n 2 - 2j 



+ 0(5n{n 2 - 2j + l)" 1 ). 



ny/j(n 2 - j) 

The estimate holds for j = n 2 /2 as well, since S 1 / 2 !!? — > oo. 



Note that l|7tj|) implies in particular the rough bound 



Pnj < 



valid for all j < n 2 . Now, to complete the proof of Theorem 3, we use this bound for 
j < agn 2 and Lemma 10 for j > a$n 2 . First 



E„ 



4(i) 



Here 



^ Pn,j + 51 p n ,j = Si + S 2 . 

j<cton 2 a.Qn 2 <i<fc 

J/ 2 \ 



si<2 £ r 1 / 2 = 0(na y2 )7 

j<Qon 2 
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and 



n 2 - 2j 



E9 < y 

The last sum is bounded above by 

1- It 



= + 0(Sn\ogn). 



a -n-2 y^l - t) 



dt = 2n^y^~a) + 0{nal /2 ). 



1 /2 

Therefore, since a 3> <51ogn, 



E n [Ay] < 2nv/a(l - a) + O(nao /2 ). 
So lf7U|) follows. Theorem 3 is proved. 



4.2 Proof of Theorem l(i) 

With our enhanced understanding of the distribution of A^ (1), we may now prove Theorem 
l(i). First we show that for individual boundary points, the tableau approaches the limit 
surface. Fix (a;, y) on the boundary of the square. By symmetry, we may assume that 
y = 0, < x < 1. Let e > 0. Denote 



a = L(x, 0) 



1-vT 



so that x = 2y/a(l — a). For any tableau T S T n , denote kr — t\nx\+i,ii an d let 
p T — kx/n 2 . We want to show that with high probability, \ftr — a\ is small. Note that kr is 
an integer representing the smallest j for which A^ > nx. Therefore nx < A^ T (1) < nx+1, 



\%r(i)-x 



< 



(77) 



The function f(t) := L(t,0) = (1 - y/l - t 2 )/2 is monotonically increasing and uniformly 
continuous on [0,1]. Choose a S > such that |£ — i'| < S implies \f(t) — f(t')\ < e/3. 
Choose numbers = clq < a\ < a2 < ■ ■ ■ < ajv = 1/2 such that a^+i — dj < e/3, 
i = 0,1,2, ...,N- 1. Denote = f^ 1 (a l ) = 2^(1 - a,). 
Let T S 7Z, be a tableau that satisfies 



-A 



|Qi" 2 J 



(1) - X, 



< 



(* = 1,2, 



(78) 



(this happens with high probability, by H67[l). Let < i < N be such that a, < (3t < cti+i- 
Then clearly 



6 . 1 A K« 2 j 

2 



, (1) < -A^(l) < -X [ T a ^ in2i (l) < x l+1 + 5 - 
n n n 2 



(79) 
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Combining this with (|77|> we get, for n > 2/8, 

xt — 8 < x < x i+ i + 8. 

Therefore 

ai - -<a = f(x) < a i+1 + -, 
and, since also < /3t < a i+1 and c^+i — a, < e/3, we get 

\[3 T - a\ < e. 

Summarizing, we have shown that 



\ (T e T n : \pr -a\<e) 

>P n (tgT„ :V* = 1,2,. 



,,7V, 



— ^ (^) — 



2 / ?l — »OG 



1. (80) 



Theorem l(i) now follows easily. It is enough to say that, because of the monotonicity 
of the tableau ti j as a function of i and j, and the monotonicity of the limit surface 



function L, given e > we can find finitely many points (xi,yi), (22,2/2) 
[0, 1] x [0, 1] \ {(0, 0), (1, 0), (0, 1), (1, 1)} such that the event inclusion 



, (xn,Vn) € 



TeT„: max 

l<i,j<n 
JV 



2 l i,J Li \ > 

71 \nn 



>e C 



2=1 l ~ 



> 



10 



(81) 



holds. But now, the P n -probability of each of the individual events in this union tends 
to as n — > 00 - because of Theorem l(ii) for the points (xi,yi) in the interior of the 
square (using the continuity of the function L), and because of (|80|l for the points on the 
boundary. ■ 



5 The hook walk and the cotransition mea- 
sure of a diagram 

In this section, we study the location of the k-th entry in the random tableau T 6 T„, 
when k w a ■ n 2 . The idea is to condition the distribution of the location of the k-th 
entry on the shape Ay of the fc-th subtableau of T. Given the shape Ay, the distribution 
of the location of the fc-th entry is exactly the so-called cotransition measure of A T (see 
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below) . We know from Theorem 8 that with high probability, the rescaled shape of A^ is 
approximately described in rotated coordinates by the level curve v — g a {u). Romik |3U| 
showed that the cotransition measure is a continuous functional on the space of continual 
Young diagrams, and derived an explicit formula for the probability density of its it- 
coordinate. By substituting the level curve g a in the formula from [3U| , we will get exactly 
the semicircle density (|2Jl, proving Theorem 2. 

Let A : A(l) > A(2) > . . . > A(m) > be a Young diagram with k = |A| = £\ X(i) 
cells. A cell c = £ A (1 < i < m, 1 < j < X(i)) is called a corner cell if removing it 
leaves a Young diagram A \ c, or in other words if j — X(i) and (i = m or X(i) > A(i + 1)). 
If T is a Young tableau of shape A, let c max (T) be the cell containing the maximal entry 
k in T. Obviously c max (T) is a corner cell of A. 

The cotransition measure of A is the probability measure fix on corner cells of A, which 
assigns to a corner cell c measure 

*<«>-^ (82) 

(with d(X) as in ©■) This is a probability measure, since one may divide up the d(X) 
tableaux of shape A according to the value of c max (T); for any corner cell c, there are pre- 
cisely d(X\c) tableaux for which c max (T) = c. In other words /i\ describes the distribution 
of c max (T), for a uniform random choice of a tableau T of shape A. 

It is fascinating that there exists a simple algorithm to sample from fi\ . This is known 
as the hook walk algorithm of Greene-Nijenhuis-Wilf, and it can be described as follows: 
Choose a cell c = £ X uniformly among all k cells. Now execute a random walk, 

replacing at each step the cell c with a new cell c', where c' is chosen uniformly among all 
cells which lie either to the right of, or (exclusive or) below c. The walk terminates when 
a corner cell is reached, and it can be shown that the probability of reaching c is given 
by (|82() . Figure 6 shows a Young diagram, its corner cells and a sample hook walk path. 

Now consider a sequence A„ : A„(l) > A„(2) > A„(3) > ... of Young diagrams for 
which, under suitable scaling, the shape converges to some limiting shape described by 
a continuous function. More precisely, let f\ n (x) be as in (J2J, and let g\ n be its rotated 
coordinate version. Let : [0,oo) — > [0, oo) be a weakly decreasing function, and let 
goo be its rotated coordinate version, a 1-Lipschitz function. In this more general setting, 
think of g\ n and g^ as functions defined on all R. Assume that: there exists an M > 
such that foo(x) = for x > M, and on [0, M] f is twice continuously differentiable, 
and its derivative is bounded away from and oo (equivalently, for some K < < K', 
3oo(w) = \u\ for u (K,K'), and g is twice continuously differentiable in \K,K'\ with 
derivative bounded awaw from -1 and 1). Finally, assume that 

\\g\ n ~ ffoolU * 0. 
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Figure 6: A Young diagram, its corners and a hook walk path 



For any n, let (I n ,J n ) be a yL\ n -distributed random vector. Let X n = I n /n,Y n = J n /n. 
We paraphrase results from [5U| . 

Theorem 10. (Romik j3U], Theorems 1(b), 6) As n — > oo, (X n ,Y n ) converges in 
distribution to the random vector 

V + U V-U \ 
2 ' 2 ) ' 

where V — goo{U) and U is a random variable on [K,K'] with density function 




with 

l-M i-K' 
A = \ f ao {x)dx = / (fl"oo(w) - \u\)du 

JO JK 

and the integral in the exponential being a principal value integral. 

Proof of Theorem 2. We may assume that < a < 1/2. The proof of Theorem 2 
now consists of an observation, a remark, and a computation. 

The observation is that the distribution of the location of the fc n -th entry in a random 
tableau T £ 7~ n is the distribution of the maximal entry in the shape A^" of the k n - 
th subtableau of T. Because by Theorem 8, this shape (suitably rescaled and rotated) 
converges in probability to g a (Theorem 2 assumes k n /n 2 —* a), we may apply Theorem 
10 and conclude that Theorem 2 is true with a density for U a given by taking — g a , 
A = a, —K = K' = ^2a{l - a) in (Si- 
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The remark is that the above is not quite true, since g a does not satisfy the assumptions 
of Theorem 10! The problem is that 



- lW Q (-VMl-a)+e) - lim&U/2a(l -«)-£) = 1, 

e\0 e\0 

so the derivative is not bounded away from -1 and 1. However, since this only happens 
near the two boundary points, going over the computations in |3l)| shows that this is not 
a problem, and the formula (|83[) is still valid in this case 4 . 

The computation is the verification that (|83[) gives the semicircle distribution @ 
under the above substitutions. We compute, using l|33|) and the identity cos(tan _1 w) = 

(l + w 2)-l/2. 

2 2 



^{x-K){K' ~x) = VM1 — a) — x 2 



exp — / du 



2 J K x — u J y 1 — a? y 2a(l — a) — x 



2 ' 



V 2 J y y/2a(l -a)-x 2 J 

(l-4a(l-a))x 2 Y 1/2 _ V 2a (l -a)- x 2 



2a(l-a)-x 2 J 2^a(l-a)^-x 2 

Multiplying the above expressions gives 

<p u (x) = —- -^2a(l-a)~x 2 , \x\ < y/2a(l-a), 

7ra(l — a) 

as claimed. 



6 Plane partitions 

We prove Theorem 4 on the limit shape of plane partitions of an integer m over annxn 
square diagram, when n e = o{m). The basic observation relating this to Young tableaux 
is that in this asymptotic regime, almost all plane partitions have distinct parts. A plane 
partition with distinct parts can be completely described by separately giving the order 

4 Alternatively, one may use the less explicit formula (8) from jSU], which is valid even 
without the assumption that g'^ is bounded away from ±1, to verify directly that J3j is 
the cotransition measure of g a . 
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structure on its parts - a square Young tableau - and an unordered list of the parts, 
which is simply a linear partition of in into n 2 distinct parts. The structure of these linear 
partitions is described by a limit shape theorem due to Vershik and Yakubovich [37], |38| . 
Combining these results will give us our proof of Theorem 4. 

We will use a result of Erdos and Lehner on partitions into a fixed number of sum- 
mands. 

Theorem 11. (Erdos-Lehner |9j) Let p(m, k) denote the number of partitions 
of m into k parts. Let q{n, k) denote the number of partitions of m into k distinct parts. 
If m and k are sequences of integers that tend to infinity in such a way that k 3 — o(m), 
then 

q{m,k) ^ ^ 

p{m, k) 

In words, if fc 3 = o(m), almost all partitions of m into k parts have no repeated parts. 
Proof. This is a combination of Corollary 4.3 and Lemma 4.4 in [§]. ■ 



For a Young diagram A, denote by p\(m) the number of plane partitions of m of shape 
A. Denote by q\(m) the number of plane partitions of m of shape A with all parts distinct. 

Lemma 11. Let A m be a sequence of Young diagrams. Let k m = |A m |. If fc 3 ^ = o(m) 
asro-xx), then 

qx m (m) > 1 

PX m (m) m^oo 

In words, if fc 3 ^ = o(m), almost all plane partitions of shape A m have no repeated parts. 

Proof. If A is a Young diagram of size k = |A|, a plane partition of m of shape A is 
described by the order structure on its parts, and the unordered set of the parts. This 
gives the equation 

q\(m) = d(X)q(m,k). 

We claim that 

Px(m)<d(X)p(m,k). (84) 

This will prove the claim, since then we will have 

q(m,k m ) < q Xm (m) < ^ 
p(m,k m ) ~ p\ m {m) ~ 

and the Lemma will follow from Theorem 11. To prove (|84|l . we define a mapping 
that assigns injectively to each plane partition 7r = (pi,j)(i,j)e_x a P an " (2^ A*)? where 
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T = (U,j)(i,j)e\ is a Young tableau of shape A, and fi : fi(l) > fi(2) > . . . > fi(k) is 
a partition of m into k parts. The mapping is defined as follows. Define a linear order 
"-<" on the cells of A, by stipulating that 



-< (i')i') Pi,3 >Pi',j' or 



Pi,j — Pi',j' and (i < i' or (i = i' and j < j')) 



Let (iijj'i) -< («2, ia) -<...-< (ik,jk) be the cells of A sorted in this linear ordering, and 
set ij,^ = / and ^(7) = i = 1, 2, . . . ,fc. 

It is easy to verify that the mapping is injective and has the required range. See Figure 
7 for an illustration. ■ 

1 2 5 8 12 

77652 m 369 10 

T = 

7 6 5 5 4 11 13 

7T = > 

7 5 2 7 
6 

/i = 7, 7, 7, 7, 6, 6, 6, 5, 5, 5, 2, 2. 
Figure 7: Illustration of the proof of Lemma 11 



Next, we recall the Vershik-Yakubovich limit shape theorem for partitions of m into 
k distinct summands, when k = o(\fm). 



Theorem 12. (Vershik-Yakubovich |37j, |38j) Let m = m n and fc = fc„ 

grow to infinity as a function of some parameter n, in such a way that k = o(y / m). Let 
A n : A„(l) > An (2) > . . . > A„(fc) be a sequence of uniform random partitions of m into k 
distinct parts. Then for any t > 0, e > 0, 



I in 
-#{1 < I < k : A„(Z) > -0 



> e 



0. 



In words, the graph of the Young diagram of a uniform random partition of m into k 
distinct parts, when k = o(y / m), will with high probability resemble the limit shape e - *. 5 



'Actually, it is more correct to say that this is the graph of the conjugate partition A'. 



50 



Proof of Theorem 4. Let ir = (pi,j)^ =1 be the random plane partition of m over 
the square diagram □„. Since n 6 — o(m) and |D„| = rt 2 , by Lemma 11 we may assume 
that 7r was chosen uniformly among all plane partitions of m of shape □„ with all parts 
distinct, since this is a set of probability close to 1 in V n ,m- Equivalently, by the remarks at 
the beginning of this section, we may assume that 7r is selected by choosing independently 
a random Young tableau T £ T n and a random partition /i : > m(2) > . . . > n(n 2 ) of 
m into n 2 distinct parts, then setting p^j = /i(tij). 

Fix < x,y < 1. Let a = L(x, y). Let i = [nx\ + 1, j = [ n ^J + 1> an d /3 = U,j/n 2 . 

We need to show that S^(x, y) — (n 2 /m)pij is with high probability very close to 
— log a. Let e > be small. From Theorem l(i), we know that with (asymptotically) high 
probability 

\p-a\<e. (85) 

Now, apply Theorem 12 for the random partition /i with t = — log(a — 2e). This gives 
that with high probability 



< e. 



\#{l<l<n 2 : > ™ (- log(a - 2e))} - (a - 2e) 

or equivalently, since > /z(2) > . . . > fJ,(n 2 ), 

4j max (l < / < n 2 : ^(Z) > ^(- log(a - 2e))| - (a - 2e) 
n I n > 

This implies in particular that 

maxjl <l <n 2 : fj,(l) > ^-(- log(a - 2e))j < n 2 (a-2e + e) = (a-e)n 2 , 
hence, since by > (a — e)" 2 , 

= < — (- log(a - 2e)). 

Apply Theorem 12 again with t = — log(a + 2e). This gives that with high probability 

\# f 1 < I < n 2 : M (/) > ™ (- log(a + 2e))} - (a + 2e) 

or equivalently 

\ max (l < I < n 2 : > ^(- log(a + 2e)) ) - {a + 2e) 

In particular this gives that 

maxjl <l<n 2 : fj,(l) > ™(- log (a + 2e))| > n 2 (a + 2e-e) = (c 



< e, 



< e. 



= (a + ejn , 
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hence, since by (JSS), t^j < (a + e)n 2 



TTl 

Pi,j = v(U,j) > -5 (- log(a + 2e)). 



We have shown that the event 

o 

n 

- log (a + 2e) < — pi,j- < - logfa - 2e) 

TO 

holds with asymptotically high probability. Since e was arbitrary the result follows. ■ 

7 Computations for the rectangular case 

The proof of Theorem 5 involves exactly the same ideas as the proof of Theorem 1, with 
some more computations, which we include here for completeness. The proof that Theorem 
6 follows from Theorem 5 is completely identical to the proof in section 6 that Theorem 4 
follows from Theorem 1. 

Fix < 8 < 1 and < a < 1. Our starting point is the rotated-coordinate formulation 
of the variational problem whose solution will yield the a-level curve of the limit surface 
Lg. The computations leading to this variational problem are obvious generalizations of 
the corresponding computations for the square case 6 = 1, and are omitted. 

Variational problem - the rectangular case. A function h : [-0V2/2, V2/2] - 
[0, oo ) is called a-admissible if h is 1-Lipschitz, and satisfies 

h(-eV2/2) = 6V2/2, (86) 

h(V2/2) = V2/2, (87) 

-V2/2 

(h(u) - \u\)du = ad. (88) 



1-6V2/2 

Find the unique o-admissible h that minimizes 

J(h) = — / / h'(s)h'(t)\og\s-t\dsdt. 

2 J -8^2/2 J-6V2/2 

To derive the minimizer, first consider the case when a is small. In that case, we make 
an assumption on the form of the minimizer similar to (|31|l . but with a non-symmetric 
interval [-(3 x (a), f3 2 {a)], where 0i € (O,0s/2/2), (3 2 £ (0,^2/2). That is, we assume that 
h' has the form 

-1, if -0V2/2 < s < -Px, 
/''(-■) is { G(-l,l), if -0i<s<fo, (89) 
+1, if f3 2 < s < y/2/2. 
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Replace the conditions (|86[) . (|87|) . I|88l) with the equivalent set of conditions 



h{-9\/2/2) = 6V2/2, 
V2/2 

h'{u)du = (1-61)72/2, 

-6V2/2 

V2/2 1 + q2 

uh'(u)du-\ = 9a. 

-0V2/2 4 



(90) 
(91) 

(92) 



(In the square case 9 = 1 we did not impose the restriction (|91() on h as we expected h 
to be even, i.e. h' to be odd, so that the condition l|91|l would be met automatically. Not 
anymore in the rectangular case!) Then the counterpart to l|32|) is: for ,s e (— ft, ft), 

h'(t) \og\s-t\dt- Xs-fi 

+ (s + 9V2/2) log \s + eV2/2\ + (s - 9V2/2) log I* - 9V2/2\ 
- (s + ft) log \s + ft| - (s ~ ft) log \s - ft| 

+ ft-ft + (l-0)V2/2 = O. (93) 

Here A, /1 are the Lagrangian multipliers dual to the constraints (|92|l and ({91(1 respectively. 
Differentiating 193J) with respect to s gives 



-dt = X + log 



s + ft 



log 



ft - s 



se(-ft,ft). (94) 



i-lhe-t' "'"°s + 6»V2/2 %/2/2-s' 

Introduce a = (ft + ft)/2, 6 = (ft — ft)/2, and substitute s = a£ + 6, i = arj + b. The 
above equation becomes 

rl h'(ar) + b) 
-i Z-v 



-d?y = A + log + log 



here 



7i 



ft - ft + 6\/2 



72 



71 + £ 72 - r 

ft - ft + V2 



(95) 
(96) 



ft+ft ' ft+ft 

and it is easy to check that 71,72 > 1. Applying Theorem 9 to (|95|l and using Lemma 8 
we obtain 

h 'te + 6 ) = 7r 2 (1 _ 1 e2) i /2 (^ + J & 71) - 72)) + (1 _^ 2)1/2 



7T(1 -e 2 ) V2 



A - log 71 + \/7i - 1 - log 72 + \hi - 1 



tan 



(1-0(72 + 1) 



tan 



'(1- 


-0(7i 


-1) 


(1- 


K)(7i 





+ c (i~er l/2 , (97) 
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c',c being arbitrary constants. As in the symmetric case, if h'(s) is to be bounded for 
s G (—fa, fa) (i-e. for £ e (—1,1)), necessarily 



c = 0, A = log 7l + J*/* - 1 + log 72 + J 7 | - 1 



(98) 



So we have 



ti(a£ + b) = - 



tan 



'(i±%-4-tan- 
(1-0(72 + 1) 



' (1-0(71-1) 
(1+0(71 + 1) 



(99) 



for which indeed \h'(a£ + b)\ < 1 holds. 

We still need to find fa and fa. Using i|89|). rewrite (|9"T]) and l|9"2|) as, respectively, 



= fa- fa 



-0i 



tti(t)dt + -{0l+pl) 



9a. 



(100) 
(101) 



Now evaluating these integrals using (|99|l , this gives the equations 



7| - 1 - 72 



7? - 1 - 7i 



2 ^ 



1-7-i +li\h' 



1 



Kfa-fa) + ^(Pl + (3 2 2 ) 



fa -Pi 



9a, 



(a = (fa + fa)/2, b — (fa — fa)/2). Excluding fa via (|96[) . we obtain two equations for 
71,72, namely 

1 



7i - #72 = 



((7 2 2 -l) 1/2 -72)-((7 1 2 -l) 1/2 -7i) 



9a = 



(1 + 9f - ( 7l - 9 l2 ) 2 (l + i 



2(71 + 72) 2 



4(7i + 7 2 ) : 



E 



l-7i +7n/7f-i 



(102) 
(103) 



It seems a minor miracle that these equations can be solved explicitly. Here is how. 
Isolating the difference of the radicals in the first equation, multiplying both sides of the 
resulting equation by the sum of radicals and cancelling the common factor 7! + 72 , we 
obtain 

I I 1 -4- ft 

;(72 ~7i)- 



7 1 2 -l+\/72^-l 



1 



(In particular, 72 > 71.) Combining this with the initial equation, we express the radicals 
as linear combinations of 71 , 72 : 



7?-l = - 



7|-1 



1 + 9 2 
1-9 2 

29 
1-9 2 



29 



7i 
7i 



1 



r72, 



r72- 



(104) 
(105) 
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Plugging these expressions for the radicals into 1|103|) . after collecting like terms, we obtain 
a quadratic equation for x — 72 /71 : 



x 2 (9(l -a) + a)- x(l - 6)(1 - 2a) - (1 - a + 9a) = 0. 



Consequently 



Jl 

7i 



9a + 1 — a 



a + 9(1 - a)' 

which, for a < 1/2, exceeds 1. Squaring both sides of (|104|l and substituting 72 = sryi, 
we solve for 71 to obtain 



7i 



a + 9(1 - a) 
2y/6a(l - a)' 



72 



6a + 1 — a 
2y/6a(l - a)' 



(106) 



Direct checking reveals that these 71 , 72 satisfy the equations (|104|l , (|l(J5fl themselves as 
long as 

a<a*:=-^—. (107) 



For a > a* , the gammas do not satisfy (|104[> . More precisely, 71,72 would have satisfied 
this equation, had we considered the negative value of y/j^—l. However, we need the 
positive value only. The source of the trouble here is that 71 = 1 for a = a* . Tellingly, 
the boundary point (— /3i, h(— reaches the corner (-6^2/2, 6V2/2) of the rotated 
rectangle at a — a*. Using l|96|l . we obtain: for a < a* , 



P 1 = y / 26a(l - a) - a(l - 9)V2/2, 
(3 2 = y/29a(l - a) + a(l - 9)^2/2. 



Using (fTUBll and ltTU5|l . we simplify Q to 



h'(t) = - tan 

7T 



(l-0)y/a(l-a)+£V6(l-2a) 



where 



t-b _ t- a(l-9)V2/2 



(108) 



(109) 



s/29a(l - a) 



Furthermore, using 



h(s)=px+ I ti (t)dt = fix + a I h'(ar] + b)dr] 
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(O, and ||STJ) we obtain 



h(s) = h{a£ + b) = fh 

(C + 7i)tan- 1 



2a 

7T 



' (1-0(71-1) 
(1+0(71 + 1) 



+ (C-72)tan 1 



' (1+0(72-1) 
(1-0(72 + 1) 



1 / • 1 ^ T 

2 ( Sm ^ + 2 



71-1 



7f - 1 + 77(71 - 1) 



Finally, combining l|98|) and (|10€>|) . we compute 



A = log 



1-a) 



log- 



(1-a) 



y/6a{l - a) y/6a{l - a) 



log 



1 - a 



(110) 



(HI) 



the same value as in the square case! 

It remains to consider the range a* < a < 1/2. In this case, it turns out that the for- 
mulas for the parameters f3\ , (3% remain the same, while the formula for the corresponding 
h! changes slightly. What is different is that now h'(s) = 1 for s e [-9V2/2, -fa] and 

h{-px) = eV2-Px. 

The latter condition means that now the boundary point [— h{—(3x)\ lies on the longer 
side of the rectangle. The starting point now is a modification of 1|94(1 . stemming from 
h!{s) = 1, rather than —1, for s 6 [— 9^/2/2, — j3i], namely 

1 h'(a V + b) HI., 1-^ 
dr] = A — log + log ■ 



71 + e 



72 -£ 



This leads to 

k 'te + & ) = 7r 2 (1 __g 2) l/2 ~ J & 71) - ^(-0 72)) + (1 _ r al/2 : 

(c' being a constant), compare with the first line in l|97|l . whence to 



A = log 72 + J 7| - 1 - log 71 + J 7l 2 - 1 , 



compare with J^HJl, and 

^+6) = ^ 



tan 



' (1-0(71-1) 
(1 + 0(71 + 1) 



tan 



' (1 + 0(72-1) 
(1-0(72 + 1) 
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compare with i|99[l . After integration, the final formula for h(s) is 



h(s) =6V2-/3! 



2a 



+ i(sm- 1 C + |) ^7i 2 -l + ^7l-l) +|(l-7i) 



(112) 



We add that, despite the difference between the two formulas for A - the one above for 
a > a* and (|98|l for a < a* - the eventual expression is still that in 11111) . The "secret" 
is that 

— = K/o- y/6(l-a)\ 
1 2y/0a(l - a) ' 



with y/a — y/ 6(1 — a) changing its sign at a*. 

It remains to prove that h is indeed a minimizer. Let a < a*. Consider s € 
[— 6\/2/2, — 0i\. Since /i'(s) = —1, we need to show that F(s,a) > 0, where F(s,a) is 
the left-hand side expression in i|93[l . The above computations show that F(s,a) = 
for s G (— 01,02)- As in the square case, for fixed s G [— 6>/2/2, 0\] let a € (a, a*) be 
defined by s = —j3\(a). Then F(s, a) — 0, and we will prove F(s, a) > if we show that 
dF(s,x)/dx < for all x £ [a, a). Since € (— 0i(x), 02(x)), F(0,x) — 0, and we use 
the latter equation to exclude the Lagrangian multiplier /i in the expression for F(s,x). 
Then, an easy computation shows that 

dF(s,x) dhM) , , ,, f 02 dh x (t), ,., dA , oN 

where ft = 0i(x), A = A(x) are given by l|108|) and i|lll|) . Let us evaluate dh x (s)/dx for 
s 6 (—01,02)- Differentiating (|94|) with respect to x we obtain 



/. 



th dh x (t)/dx dX _ 1 



s — t dx x(l — x) 

Then, using Theorem 9 and (|36() . 

dh x (s) £ c , . 

— — — = 1 (114) 

<9x 7rx(l-x)(l-e 2 ) 1/2 (l-e 2 ) 1/2 ' 

Here we have to set c = 0, as the equation 1100(1 - upon differentiation with respect to x 



leads to 



Hence 

dh x (s) 



1 1 ^ = °- < 118 ) 

(116) 



dx 7rx(l-x)(l-^ 2 ) 1 /2- 
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Plugging this expression into (|113f) . integrating by parts, and using 14U|) . we transform 
itn^I into 

dx x(l-x) v ; 

Let a E (a*, 1/2], and s < —/3i(a) again. Since now h'(s) = 1, we need to show that 
F(s,a) < 0. Let a S (a*, a) be defined by — s — fli(a). (a exists, uniquely, because Pi(x) 
is decreasing on (a, a*) and — s < (3i(a).) Then F(s,a) = 0, and so again we need to 
show that dF(s,x) / dx < 0. The formula 1114(1 continues to hold, and so does (|115fl . since 
now we have 

0% 

h'(t)dt = P2 + P1-6V2/2, 

-Pi 

and h'(—(3i) — 1. Therefore (|116fl remains valid, which implies (|117|l . 

Analogously, F(s,a) < for s > ^(a) and a £ (0, 1/2]. This finishes the proof that 
h is the minimizcr, the claim which forms the core of the proof of Theorem 5. ■ 



8 Discussion 

8.1 Plancherel measure 

Let 34 denote the set of Young diagrams of area k. The Plancherel measures are the 
family of probability measures /it on 34, defined by 

MA) = ^T' ( Ae ^)- ( 118 ) 

Alternatively, fik is sometimes defined as a measure on all Young tableaux of size fc, where 

d(shape(T)) 



k\ 



The measure on diagrams is then the projection of the measure on tableaux under the 
mapping that assigns to each tableau its shape. The measures [ik are a 'projective family of 
measures, in the following sense: If T is a /ifc-random tableau, then the tableau T' of size 
k — 1 obtained by deleting the A:-th entry from T is a /Ufc_i-random tableau. Therefore, all 
the /Life's can be encompassed by a single object P, the infinite Plancherel measure, which 
is a measure on infinite tableaux - i.e. fillings of the squares in the positive quadrant 
of the plane with the positive integers that are increasing along rows and columns - for 
which the marginal distribution of the shape of the k-th subtableau (the set of squares 
where the entry of the infinite tableau is < k) is given by 1)118(1 . In other words, P can 
be thought of as a measure on all sequences = Ao C Ai C A2 C ... of Young diagrams, 



58 



where Xk has size k and is obtained from Xk-i by the addition of a box. So P is simply 
a natural Markovian coupling of the measures (|118|l . known sometimes as the Plancherel 
growth process. 

Much is known about Plancherel measure. It arises naturally in representation theory, 
as a natural measure on the irreducible representations of the symmetric group, and in 
combinatorics, as the distribution of the output of the RSK algorithm applied to a uniform 
random permutation in Sk- In particular, the length of the first row Afe(l) of a /i^-random 
Young diagram has the same distribution as the length l n (Tr) of the longest increasing 
subsequence of a uniform random permutation tt in Sk, an important permutation statistic. 

Logan-Shepp (20j and Vershik-Kerov jSS], EH] proved that the graph of a /Zfe-random 
Young diagram, when rescaled by a factor of \fk along each axis and drawn in rotated 
coordinates, with high probability resembles the limit shape 



f (usmT 1 (u/V2) + V2 - u 2 ) M<V2, 
\u\ \u\ > V2, 



see Figure 8. 




-1.5 -1 



-0.5 



0.5 1 



1 . 5 



Figure 8: The limit shape v = Q(u) 



Gribov ^3] noted that this can be reinterpreted as a theorem on the limit surface of 
the Plancherel-random tableau, much in the same spirit as Theorem 1. If T is a /^-random 
Young tableau, then after rescaling the graph of T is approximately described in rotated 
coordinates by the surface S : D — > [0, 1], where 

D = {{u,v) : \u\ <V2, \u\ < v < fl{u)} 

is the two dimensional domain bounded between the graphs |u| and and for each 
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< a < 1, the a-level curve of £ is 

{(u,v) G £> : |u| < V2^, u = Vafi(!i/Va)}, 



a shrunken copy of 17. 

The approach in the papers of Logan-Shepp and Vershik-Kerov was the variational 
approach, of analyzing the limiting integral functional arising from l|118|l . Kerov j!8j 
considered the following more dynamical approach: Assume that we have selected the 
//fc-random diagram Afc. Since under P, the sequence Ai C A2 C ... is a (nonhomogeneous) 
Markov chain with values in the Young graph, there is a measure v on the exterior corners 
of Afc (the boxes in the complement of Afc that can be added to Afc to form a Young diagram 
of size k + 1), such that if we choose a z/-random corner of Afc and add the new box there, 
the resulting diagram Afc + i will have distribution fik+i- lu other words, v is the probability 
transition measure of the Markov chain (Afc). It is known as the transition measure of the 
diagram A^, and is in a sense dual to the co-transition measure discussed in section 5. 

Kerov showed that in the limit when the graph of the diagram Afc becomes a smooth 
curve, the transition measure converges to a limit. Imagine that in the limit, instead of 
attaching a new box at a z^-random corner, one attaches a ^-fraction of a box at each 
corner. So the curve grows in the "tangent" direction given by v. Thus, the Plancherel 
growth process can be described in the limit as a smooth flow on the (infinite-dimensional) 
space of shapes. Kerov showed that £l(u) is the unique shape which, after rescaling, is 
invariant under this flow, and that this fixed point is an attractor of the flow; this explains, 
in a way, (though does not formally prove) its appearance as the limit shape for Plancherel- 
random diagrams. Remarkably, the transition measure of £1 (the limiting direction of the 
flow) is the semicircle distribution. 

Another interesting direction stemming from the study of Plancherel measure is the 
connection to longest increasing subsequences of random permutations. The limit shape 
result of Logan-Shepp and Vershik-Kerov implies that the length Z„(7r) of the longest 
increasing subsequence of a random permutation n 6 S n is with high probability at least 
(1 — o(l))2y / n. Using additional arguments (which were an inspiration for our proof of 
Theorem 3), Vershik and Kerov showed also that l n (ir) is with high probability at most 
(1 + o(l))2y / n, solving the so-called Ulam's problem. More recently, Baik, Deift and 
Johansson 3 showed that the fluctuation of l n (ir) around its asymptotic value has a 
limiting distribution. More precisely, 



in distribution 



Here F is the Tracy- Widom distribution from random matrix theory, denned as 
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where u(x) is the solution of the Painleve II equation u"{x) — 2u(x) 3 + xu{x) that is 
asymptotic to the Airy function Ai(x) as x — ■» oo. Other results along those lines can be 
found in [B|. [T7|. see also the survey [T] 

The distribution F appears in random matrix theory as the limiting distribution of 
the maximal eigenvalue of a GUE random matrix. Following the Baik-Deift- Johansson 
result, it was found that there are many striking parallels between Plancherel measure 
and random matrix ensembles, see jS], |23| . |16|. In particular, the transition measure of 
the Plancherel-random diagram converges to the semicircle law, which is also the limiting 
distribution of the empirical eigenvalue distribution in the GUE and GOE random matrix 
ensembles. Ivanov and Olshanski ^H] showed that this similarity is no mere coincidence, 
but in fact appears also in the finer fluctuations of the transition measure and eigenvalue 
distribution measure around the semicircle distribution. 

8.2 The random square tableau as a deformation of 
Plancherel measure 

The reader familiar with the works of Logan-Shepp and Vershik-Kerov will undoubtedly 
have noticed the similarity between these results and our analysis of the square tableau 
model. Define for each positive integer n and each 1 < k < n 2 , the probability measure 
v n ,k on y k , by 

^<A> = «™ ( 119 ) 

where d(O n \ A) is taken as if A ^ The measure f n ,fc 1S the distribution of the fc-th 
subtableau of a random n x n square tableau, and our entire approach revolved around 
the analysis of its properties. It is remarkable how many of the ideas used in the study 
of Plancherel measure we have found useful in our study of square tableaux; first, and 
most obviously, the variational problem that arises from l|119fl resembles the variational 
problem studied by Logan-Shepp and Vershik-Kerov. Although our approach in solving 
the variational problem relied on the more methodical use of the inversion formula for 
Hilbert transforms (an approach that could be applied the Plancherel case as well!), we 
were greatly inspired by the methods used in the Plancherel case. Secondly, our proof of 
Theorem 3 and the treatment of the boundary of the square also follows closely the ideas of 
Vershik and Kerov (with the notable difference, that our proof of the upper bound uses the 
lower bound!). Finally, our Theorem 2 on the location of particular entries, was inspired 
by Kerov's differential model JS] for Plancherel growth. By postulating the existence of 
an analogous differential growth model for the fc-subtableaux of the square tableau, we 
were able to guess Theorem 2 from the formula (|66ll . This was later verified by a different 
method, using the result from |3(J) . 
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Take another look at (|119f) and H118fl . The defining equations for fik and is n ,k seem 
superficially similar at best. In fact, they are closely related, and when k is very small 
these measures are quite close. To make this precise, we first note the following curious 
identity. Define the falling power ar of a as a^ b = a(a — 1) (a — 2) ... (a — b + 1). Then: 

Lemma 12. If A e y k , A c □„, then 

VnAX) _ (n + j-l)U. n?iV (n + J- 1)^ 

MA) (n 2 )i fe 

Proof. Use the hook formula 10 • A computation similar to the one in the proof of 
(t?2*jl shows that many of the terms cancel. We omit the relatively simple details. ■ 

It follows using elementary estimates, which we again omit for the sake of brevity, that 
Theorem 13. If n — > oo and k = k{n) is such that k = o(n 2 / 3 ), then 

^n,fc(A) ^ 

uniformly on the support y n ,k of v n ,k (the set of diagrams of size k contained in □„). In 
particular, the total variation distance 

IK.fc - Mfelli := } y l^n,fc(A) - MfcWI * °- ■ 

^ — ^ n — ►oo ■ 

So in fact, when fc is small, v n ^ is a kind of deformation of the Plancherel measure 
/ifc. In particular, for k fixed and n going to infinity this implies the not-entirely-trivial 
fact that fj,k is a probability measure. We remark that other deformations of Plancherel 
measure have been used as a means to study Plancherel measure itself - see, e.g., |17j . 
The phenomenon that a small subtableau of a large random tableau has approximately 
the Plancherel distribution was observed also in [221 ( see a ^ so for related results) for 
a random tableau chosen uniformly among all tableaux of size k. Recently, it was shown 
26 that the footprint of the k tallest stacks in a random unrestricted plane partition of 
high volume also has in the limit the Plancherel distribution. 

Another related observation is the easily checked fact that 

\/a(l-a) ■ g a = —— * fi(tt), 

y ^Ja{l - a) J <*\0 

i.e. the shape of the level curves of our limit surface L converges after rescaling to the 
Plancherel limit curve f2, as one approaches the corner of the square. This is consistent 
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with Theorem 13, although is not formally implied by it, as here k is a small constant 
times n 2 . It seems likely that in the regime when k grows like o(n 2 ), but much faster than 
n 2 ^ 3 , v n ,k and fi^ become mutually singular, even though the limit shapes coincide. 



8.3 The probability of a square plane partition to have 
all parts distinct 

Denote by M n jq the total number of n x n square plane partitions of N. Let M n N be the 
total number of those partitions with all parts distinct. From Lemma 11 it follows that if 

n 6 

lim — = 0, (120) 

then 

lim ^ = 1. (121) 

n,iV-K» M nN 

Our goal is to show that (|120f> is essentially necessary for (|121|l . To motivate the statement, 
notice that the k-th largest part in a partition of N into n 2 distinct parts is n — k + 1, at 
least. So M nN = unless N > n 2 (n 2 + l)/2. 

Theorem 14. Suppose that n 4 /N -> 0. (i) If limn 6 /TV = oo, then 

lim 4^=0. (122) 

(ii) If limn 6 /N = a G (0, oo) then 

lim ^=e"«/ 4 . (123) 

n,N— >oo M n N 

Note. Thus the reduction to the plane partitions with distinct parts used in the proof 
of Theorem 4 is valid if and only if ri 6 = o(m). 

Proof sketch of Theorem 14. We prove (|122l) . 1123|) by determining the asymp- 
totic expressions of M n N and M„n- 

Part 1. Begin with A4 n N- As in the proof of Lemma 10, we notice that - given a 
linear partition of N into n 2 distinct parts - the number of the n x n square (descending) 
arrangements of these parts equals the total number ofnxn square Young tableaux. So, 
denoting by p n jy the total number of all such linear partitions, and by d(U n ) the number 
of all such tableaux, we obtain 

M nN = p nN ■ d(n n ). (124) 
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Using the hook formula Euler's summation formula for X)s=i( n — s ) ^°s( n — s )j an( i 
two identities for the Gamma function (see Bateman [5], Section 1.9) 



n „ n 

^ slogs = / logr(a;)da; + 



logT(x) 



1 



n(n +1) n 



1 



x ~ 77 logx-2; + -log27r 



+ 



log2?r, 



we obtain 



d(D n ) ~ n 11/12 V2^exp ( n 2 logn + n 2 (-21og2 + 1/2) - \ + - C 
\ 6 12 



where 



C := 



1 t 

2 ~ 12 



i~ 2 e~* di. 



(125) 
(126) 



(A cruder formula 

<f(D„) - exp(n 2 logn + n 2 (-2 log 2 + 1/2) + O(rclogrc)) 

was obtained, implicitly, in the proof of Lemma 1.) 

As for p n N, the total number of partitions of ./V into n 2 distinct parts, it is given by 



PnN = lq N t n2 ]l[(l + q e t). 



(127) 



From a more general theorem of Vershik and Yakubovich based on l|127jl . Fristedt's 
conditioning defice ^21j and an attendant local limit theorem result, it follows that 



PnN 



2ttN 



exp n log 



N 



2n z 



4N 



(1 + 0(n 4 /N)) 



(128) 



Combining l|125|) and Q 128)1 . we arrive at 



Mr 



,11/12 



N 



exp(n 2 log{N/n 3 ) + n 2 (-2 log 2 + 5/2) 

4 , 



c* 



2ttN 

~(n 6 /AN)(l + 0(n 4 /N)) + C*), 
1 log2 



6 12 
with C defined in (T2"6) . 



(129) 
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Part 2. Turn now to M n N, the number of all square n x n plane partitions of N. By 
the MacMahon formula for the number of plane partitions with at most n rows and n 
columns, 

oo 

M nN = [q N - 2 ] f[(l q e y e Y[(l qY l - n) U (1 - ^) 2 "^ (130) 

fcl i>n l>2n 

(Alternatively, this formula follows from the hook expression for the generating functions of 
plane partitions with a given shape discovered by Stanley J2H].) We will use the techniques 
from [2HI inspired by Freiman's derivation of the main part of Hardy- Ramanuj an formula 
for the (linear) partition function, see Postnikov |29j . 

Let us take a close look at the generating function in i|13U|) . which we denote p n (q). 
Set q = e~ u , Re u > 0. Taking logarithms, using 

log(l - e-™)- 1 =E^— > ( 131 ) 



;>i 3 



and changing the summation order, we obtain 



_ , i P «j 
logp„(e-") = u £ - (—3^(1 ~ e" UnJ ) 2 

e~ uru 3 ij)(unj) 



oo 

2 



E— — + un3 E 



£ (1 - e-^f +uJ2Huj)(l - er^f- 

3=1 3 3=1 

e z 1 e~ z 

z(e z -l) 2 ~ ^ + vTz 

-z\2 „2 -2 



i>{z) := (\-e- z Y -z z e~ z . (132) 

Using (|131|l . read backward, for the first and the third sums, and Euler's summation 
formula, with m = 1, for the second and the fourth sum, we obtain from (|132|l : if |u| 2 n 4 — > 
0, then 

n( -u, exp(n 2 (3/2-21 g2) + (log2)/12 + J) „ a> 
Pn\e ) ~ ^ 1/12 ■ ^ e J 

p OO 

D ;= cj)(x)dx. (133) 
j 

(According to Maple the integral of ^>(a;) equals 3/2 — 2 log 2.) Now, by l|13U|) . and the 
Cauchy formula 

MnN = [^"" 2 K(5) = ^ / ^ efe, 
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where p € (0, 1). In light of the last two equations, we set p = e u ° , and select uq that 
minimizes — n 2 log(l — er nu ) + (N — n 2 )u, that is 

u = n^ 1 log I 1 + 



N -n 2 J N ' 

Clearly u 2 n 4 -> 0. Consider |0| < 9 n = n~ 2 e„, e„ = (n 4 /iV) 1/2 , so that |6f?i 4 -> as well. 
It can be shown without much difficulty that 

1 /■ 1 f e"( N ~ n ^ , 

Gw = - — : / -, rr «w 



27TI J (l_ e -™)n 2 z lV-n 2 +l 27TI J (1 ~ e~ nv ) 

\8\<e n u -i9„ 

= A- n (1 + 0(n /N)) ~ —7= — ( — r- ) • (134) 
n" 2 (n 2 -l)P v 7 V2^N\n3J v ' 

(The last integral, extended to the closed contour obtained by connecting the points uo ± 
with a circular arc centered at the origin, is exactly 



J_[^-i] (1 _ nt) -(N- n > + n + 2 )/n = n -n* (- n " {N ~ f + n + 2) ) (_„)r 



and the supplementary integral is less than this quantity by a factor (uo/9 n ) n ~ {n A /N) n I 2 .) 
So, using H331, JT3U), 

27TI 



96(-9ra>«n 



11/12 / jy 

exp n 2 log — + n 2 (5/2 - 2 log 2) + (log 2)/12 + D + o(l) ) . (135) 



2tt7V 

The proof that the contribution of ^ [— n ,9 n ] is negligible compared with the last 
expression is based on cruder estimates, not unlike those in |26j . and we omit it. Therefore 

?1 n/i2 / N \ 
M nN ~ _ exp n 2 log — + ?i 2 (5/2 - 2 log 2) + (log 2)/12 + D . (136) 
V27r7V V n' 3 J 



Comparing (|136|) and (|129J) . we see that 

= exp(- (n 6 /4N){l + 0(n 4 /N)) + A), (137) 



M nN 



M, 



N 
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where, recalling the definition of C and D, 
A = -l-C-D 



Thus we have 



1 1 1, , , 

- - + r t e- 



e* - 1 t 2 



1 



i(e* - l) 2 i 3 



1 e-*\ , 

— -7T + — -rr ) dt 



dt V <(e* - 1) 2t 2 2t 2 



lim — 



1 e 



0. 



which proves Theorem 14(i),(ii 



6 no V *(e* - 1) 2t 2 2t 2 



= exp(- (n 6 /4^)(l + 0(n 4 /A^))), 



8.4 Open problems 

We conclude with some open problems. 

• Gaussian fluctuations. Prove a central limit theorem for the fluctuations of 
<7^ L „„2j around g a , and for the fluctuations of the cotransition measure of A^*™ 
around the semicircle distribution, in the spirit of |16j . 

• Limiting distribution of l n ,k{^)- Find a scaling sequence a n and a distribution 
function F such that, in the notation of Theorem 3, 

1"n,\_ctn 2 \ 2-^/^(1 Ql)ti in distribution ^ 
f ' 

• Limit surface for random Young tableaux of given shape. Prove a limit 
surface theorem for random Young tableaux of other shapes. In general, one can 
consider any decreasing function / : [0, oo) — * [0, oo) such that f(x)dx = 1 as 
a continual Young diagram, i.e. as a limit of the rescaled graphs of a sequence of 
Young diagrams of increasing sizes. We conjecture that for each such continual 
diagram /, there should exist a limit surface Lf, defined on the domain 

D f :={ix,y):x>0, 0<y<f(x)} 

bounded between the z-axis and the graph of /, that describes the asymptotic 
behavior of almost all random Young tableaux of shape approximated by /. 
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